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' Summary. The lectures are centered around three selected topics ol quantum 

^ ^ ' chaos: the Selberg trace formula, the two-point spectral correlation functions of 

Riemann zeta function zeros, and of the Laplace-Beltrami operator for the modular 
' group. The lectures cover a wide range of quantum chaos applications and can serve 
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as a non-formal introduction to mathematical methods of quantum chaos. 



Introduction 



Quantum chaos is a nickname for the investigation of quantum systems which 
do not permit exact solutions. The absence of explicit formulas means that 
underlying problems are so complicated that they cannot be expressed in 
\^ , terms of known (~ simple) functions. The class of non-soluble systems is very 

' large and practically any model (except a small set of completely integrable 

, systems) belongs to it. An extreme case of quantum non-soluble problems 

appears naturally when one considers the quantization of classically chaotic 
systems which explains the word 'chaos' in the title. 
' As, by definition, for complex systems exact solutions are not possible, new 

H . analytical approaches were developed within the quantum chaos. First, one 

may find relations between different non-integrable models, hoping that for 
certain questions a problem will be more tractable than another. Second, one 
^ ^ ' considers, instead of exact quantities, the calculation of their smoothed values. 

, In many cases such coarse graining appears naturally in experimental settings 

^ ' and, usually, it is more easy to treat. Third, one tries to understand statistical 

5^ , properties of quantum quantities by organizing them in suitable ensembles. 

An advantage of such approach is that many different models may statistically 
be indistinguishable which leads to the notion of statistical universality. 

The ideas and methods of quantum chaos are not restricted only to quan- 
tum models. They can equally well be applied to any problem whose analytical 
solution either is not possible or is very complicated. One of the most spec- 
tacular examples of such interrelations is the application of quantum chaos 
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to number theory, in particular, to the zeros of the Riemann zeta function. 
Though a hypothetical quantum-like system whose eigenvalues coincide with 
the imaginary part of Riemann zeta function zeros is not (yet!) found, the Rie- 
mann zeta function is, in many aspects, similar to dynamical zeta functions 
and the investigation of such relations already mutually enriched both quan- 
tum chaos and number theory (see e.g. the calculation by Keating and Snaith 
moments of the Riemann zeta function using the random matrix theory 43 ). 

The topics of these lectures were chosen specially to emphasize the inter- 
play between physics and mathematics which is typical for quantum chaos. 

In Chap.Hdifferent types of trace formulas are discussed. The main atten- 
tion is given to the derivation of the Selberg trace formula which relates the 
spectral density of automorphic Laplacian on hyperbolic surfaces generated 
by discrete groups with classical periodic orbits for the free motion on these 
surfaces. This question is rarely discussed in the physical literature but is of 
general interest because it is the only case where the trace formula is exact 
and not only a leading semiclassical contribution as for general dynamical 
systems. Short derivations of trace formulas for dynamical systems and for 
the Riemann zeta function zeros are also presented in this Chapter. 

According to the well-known conjecture statistical properties of eigen- 
values of energies of quantum chaotic systems are described by standard ran- 
dom matrix ensembles depending only on system symmetries. In Chap. |2] we 
discuss analytical methods of confirmation of this conjecture. The largest part 
of this Chapter is devoted to a heuristic derivation of the 'exact' two-point 
correlation function for the Riemann zeros. The derivation is based on the 
Hardy-Littlewood conjecture about the distribution of prime pairs which is 
also reviewed. The resulting formula agrees very well with numerical calcula- 
tions of Odlyzko. 

In Chap. 01 a special class of dynamical systems is considered, namely, 
hyperbolic surfaces generated by arithmetic groups. Though from viewpoint 
of classical mechanics these models are the best known examples of classical 
chaos, their spectral statistics are close to the Poisson statistics typical for 
integrable models. The reason for this unexpected behavior is found to be 
related with exponential degeneracies of periodic orbit lengths characteristic 
for arithmetical systems. The case of the modular group is considered in de- 
tails and the exact expression for the two-point correlation function for this 
problem is derived. 

To be accessible for physics students the lectures are written in a non- 
formal manner. In many cases analogies are used instead of theorems and 
complicated mathematical notions are illustrated by simple examples. 
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Trace Formulas 



Different types of trace formula are the cornerstone of quantum cfiaos. Trace 
formulas relate quantum properties of a system with its classical counterparts. 
In the simplest and widely used case the trace formula expresses the quantum 
density of states through a sum over periodic orbits and each term in this sum 
can be calculated from pure classical mechanics. 

In general, dynamical trace formulas represent only the leading term of 
the semiclassical expansion in powers of h. The computation of other terms is 
possible though quite tedious pp. The noticeable exception is the free motion 
on constant negative curvature surfaces generated by discrete groups where 
the trace formula (called the Selberg trace formula) is exact. The derivation 
of this formula is the main goal of this Section. 

For clarity, in Sect. Q the simplest case of the rectangular billiard is briefly 
considered and the trace formula for this system is derived. The derivation 
is presented in a manner which permits to generalize it to the Selberg case 
of constant negative curvature surfaces generated by discrete groups which 
is considered in details in Sect. El In Sects. 01 and 01 the derivations of the 
trace formula for, respectively, classically integrable and chaotic systems are 
presented. In Sect. [3 it is demonstrated that the density of Riemann zeta 
function zeros can be written as a sort of trace formula where the role of 
periodic orbits is played by prime numbers. Section |B1 is the summary of this 
Chapter. 

1 Plane Rectangular Billiard 

To clarify the derivation of trace formulas let us consider in details a very 
simple example, namely, the computation of the energy spectrum for the plane 
rectangular billiard with periodic boundary conditions. 
This problem consists of solving the equation 
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where A — /dx^ + /dy^ is the usual two-dimensional Laplacian with 
periodic boundary conditions 

^^{x + a,y)^^n(x,y + b)^^^{x,y) (2) 

where a and 6 are sizes of the rectangle. 
The plane wave 

is an admissible solution of Boundary conditions jSJ determine the allowed 
values of the momentum k 

ki = — ni, K2 — —n2 , 
a 

with ni,n2 — 0,±1,±2, . . ., and, consequently, energy eigenvalues are 

27r \^ /27r 



The first step of construction of trace formulas is to consider instead of indi- 
vidual eigenvalues their density defined as the sum over all eigenvalues which 
explains the word 'trace' 

d{E)^ ^ <5(i? (4) 

ni ,n2— — C30 

To transforms this and similar expressions into a convenient form one often 
uses the Poisson summation formula 

E /N= E / e2-'""/(n)dn . (5) 



n— — oo m— — oo 



An informal proof of this identity can, for example, be done as follows. 
First 

+00 /• + 00 

E = / fi^)9{x)^^ 

n— — OQ — oo 

where g{x) is the periodic (5- function 

+00 

9{x)= E ^(a;-"-)- 



As any periodic function with period 1, g{x) can be expanded into the Fourier 
series 

g{x)^ e^-^-c^. 

m——oo 
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CoeiEcients are obtained by the integration of g{x) over one period 

'1/2 
1/2 



Cm^ /^'^%(y)e-2™My = 1 

J-l/2 



which gives ©. 

By applying the Poisson summation formula (|3J) to the density of states 
(0)) one gets 



d{E) = 



nil ,m2 — — oo 



g27ri(mini+m2n2) ^ 



2tt V /27r ^ 



X 6 \ E - [ —nij - (^-^"2 ) I dnidn2 . 

Perform the following substitutions: E — k'^, ni — arcosip/2Tr, and n2 
brsmip/2Tr. Then dnidn2 — abrdrdip / {2Tr)'^ and 

d{E) = ^l^ ^ / /"e'(™i"™"^+'"^''™^)''(5(fc2-r2)rdrdv? 

■mi,m2— — C30 

2(27r)2 ^ L 

^ ' .m2 — — OO ^ 

mi, 7712— —oo 

where IJ,{D) — ah is the area of the rectangle, 

2 i'27r 



is the Bessel function of order zero (see e.g. 13^, Vol. 2, Sect. 7), and 



Lp = y/ {miaY + (m26)2 

is (as it is easy to check) the length of a periodic orbit in the rectangle with 
periodic boundary conditions. 

Separating the term with mi = TO2 — one concludes that the eigenvalue 
density of the rectangle with periodic boundary conditions can be written as 
the sum of two terms 

d{E) = d{E) + d^°"'\E) , 

where 

m ^ 4^ (6) 
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is the smooth part of the density and 

rf^°^='(i?) = 4^E'^o(fcLp), (7) 

^ P.O. 

is the osciUating part equal to a sum over aU periodic orbits in the rectangle. 
As 

Jq{z) J — cos (z ~ ^ ) 

V TTZ V 4/ 

the oscillating part of the level density in the semiclassical limit k ^ oo takes 
the form 

Sosc)^E) = 4^ y cos (kL, - ^) . (8) 

Let us repeat the main steps which lead to this trace formula. One starts 
with an explicit formula (like which expresses eigenvalues as a function 
of integers. Using the Poisson summation formula Q the density of states (0} 
is transformed into a sum over periodic orbits. In Sect.|21it will be demon- 
strated that exactly this method can be applied for any integrable systems in 
the semiclassical limit where eigenvalues can be approximated by the WKB 
formulas. 

More Refined Approach 

The above method of deriving the trace formula for the rectangular billiard 
can be applied only if one knows an explicit expression for eigenvalues. For 
chaotic systems this is not possible and another method has to be used. 
Assume that one has to solve the equation 

(£;„ - ff)tf'„(x) = 

for a certain problem with a Hamiltonian H . Under quite general conditions 
eigenfunctions tf'n(x) can be chosen orthogonal 



J lfn(x)!f^(x)dx S„ 



and they form a complete system of functions 

5]tf'„(x).Z^*(y) = 5(x-y). 

n 

The Green function of the problem, by definition, obeys the equation 

and the same boundary conditions as the original eigenfunctions. Its explicit 
form can formally be written through exact eigenfunctions and eigenvalues as 
follows 

The -|-ie prescription determines the so-called retarded Green function. 
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Example 

To get used to Green functions let us consider in details the calculation of the 
Green function for the free motion in /-dimensional Euclidean space. This 
Green function obeys the free equation 

{E + h^A)GP{x,y)^5{^^y). (10) 

Let us look for the solution of the above equation in the form G^-* (x, y) = G{r) 
where r = |x — y| is the distance between two points. 

Simple calculations shows that for r ^ G{r) obeys the equation 

d^G , /-IdG , 



dr^ r dr h 



2 



where E = k^. 

After the substitution 



G(r) = ri-^/2g ( 



one gets for g{z) the Bessel equation (see e.g. [221, Vol. 2, Sect. 7) 

dz^ z dz \ z 



71^+ (1--). 9 = (11) 



with = I//2- 1|. 

There are many solutions of this equation. The above -l-ie prescription 
means that when k ^ k + ie with a positive e the Green function has to 
decrease at large distances. It is easy to see that G(r) is proportional to 
g±ifcr/fi g^-j- ]^g^j.gg ^_ ^jj^g prescription selects a solution which behaves at 
infinity like 6+''^'"/'' with positive k. The required solution of I|1H) is the first 
Hankel function (see [SI, Vol. 2, Sect. 7) 

g(z) = G;ilW(z) (12) 

where G/ is a constant and H^\z) has the following asymptotics for large 
and small z 

\ TTZ 

and 

^ > \ 2ilnz/7r , ^^2 ■ 

The overall factor in p2(l has to be computed from the requirement that the 
Green function will give the correct 5-function contribution in the right hand 
side of H10|l . This term can appear only in the result of differentiation of the 
Green function at small r where it has the following behaviour 
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G(r)^Go(r)=A^r2-/ 

with 

^ ^ ink'' 

One should have 

h^AGoir) = Sir) . (13) 

Muhiplying this equality by a suitable test function /(r) quickly decreasing 
at infinity one has 

J f{r)AG„{r)dr = f{0) . 
Integrating by parts one obtains 

/ 7^/W7^Go(r)dr^-/(0). 

As both functions /(r) and Go{r) depend only on the modulus of r one finally 
finds 

where 5'/_i is the volume of the (/ — l)-dimensional sphere xl + . . . + x'j = 1. 
Using (|13|l one concludes that in order to give the (5- function term Af has to 
obey 

h'Af{f-2)Sf-^ = -l. 
One of the simplest method of calculation of Sf-i is the following identity 

/oo poo 
e"^2dx2 . . . / e^^/da;/ = tt^^^ . 
-oo J —OQ 

By changing Cartesian coordinates in the left hand side to hyper-spherical 
ones we obtain 



/>oo 

Jo 



which gives 

r(//2) 

where r{x) is the usual gamma-function (see e.g. [35], Vol. 1, Sect. 1). 

Combining together all terms and using the relation xr{x) = r(x + 1) one 
gets the explicit expression for the free Green function in / dimensions 

where v — |//2 — 1|. In particular, in the two-dimensional Euclidean space 
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GL"^(x,y) = ^i/«Q|x-y|) . (15) 

Another method of calculation of the free Green function is based on Q which 
for the free motion is equivalent to the Fourier expansion 

Performing angular integration one obtains the same formulas as above. 



The knowledge of the Green function permits to calculate practically all 
quantum mechanical quantities. In particular, using 

Im ~ttS(x) 

x + is ^ ' 

one gets that the eigenvalue density is expressed through the exact Green 
function as follows 

d{E) = --Im ( G'£;(x,x)dx . (17) 

Jd 

This general expression is the starting point of all trace formulas. 

For the above model of the rectangle with periodic boundary conditions 
the exact Green function has to obey 

+ ^ + E)Ge{x, y- x\ y') = 5{x - x')5{y - y') (18) 

and the periodic boundary conditions 

Ge{x + na,y + mb; x' , y') = Ge{x, y; x\ y) (19) 

for all integer m and n. 

The fact important for us later is that the rectangular billiard with periodic 
boundary conditions can be considered as the result of the factorization of the 
whole plane (x, y) with respect to the group of integer translations 

a; ^ X + na, y y + mb (20) 

with integer m and n. 

The factorization of the plan (x, y) with respect to these transformations 
means two things. First, any two points connected by a group transformation 
is considered as one point. Hence p9(l fulfilled. Second, inside the rectangle 
there is no points which are connected by these transformations. In mathe- 
matical language the rectangle with sizes (a, b) is the fundamental domain of 
the group 12011 ■ 
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Correspondingly, the exact Green function for the rectangular billiard with 
periodic boundary conditions equals the sum of the free Green function over 
all elements of the group of integer translations H2U|I 

C30 

GE{x,y]x\y)^ ^ G%\x + na,y + mb]x ,y') . 

n,m— — oo 

Here G^'' (x, x') is the Green function corresponding to the free motion with- 
out periodic boundary conditions. To prove formally that it is really the exact 
Green function one has to note that (i) it obeys H18|l because each term in the 
sum obeys it, (ii) it obeys boundary conditions (|19() by construction (provided 
the sum converges), and (iii) inside the initial rectangle only identity term can 
produce a i5-function contribution required in (|18(l because all other terms will 
give (^-functions outside the rectangle. 

The next steps are straightforward. The free Green function for the two- 
dimensional Euclidean plane has the form H15|l . From H17|l it follows that the 
eigenvalue density for the rectangular billiard is 

d{E)^~-lm I G£;(x,x)dx = 
71" Jd 

Im i/^^^ (fcV(ma)2 + (n6)2^ dx = 



47r ^ 



D 



= 4^ + ^E' MkL,) (21) 

which coincides exactly with JBJ and |7J obtained directly from the knowledge 
of the eigenvalues. 

The principal drawback of all trace formulas is that the sum over periodic 
orbits does not converge. Even the sum of the squares diverges. The simplest 
way to treat this problem is to multiply both sides of H21|l by a suitable test 
function h{E) and integrate them over E. In this manner one obtains 



E 



h{E^) = l^ r hiE)dE + ^Yl r h{E)MVELp)dE 
47r Jo 47r ^Jq 



When the Fourier harmonics of h{E) decrease quickly the sum over periodic 
orbits converges and this expression constitutes a mathematically well de- 
fined trace formula. Nevertheless for approximate calculations of eigenvalues 
of energies one can still use 'naive' trace formulas by introducing a cut-off on 
periodic orbit sum. For example, in Fig. ^ the result of numerical application 
of the above trace formula is presented. In performing this calculation one 
uses the asymptotic form of the oscillating part of the density of state (jHJ 
with only 250 first periodic orbits. Though additional oscillations are clearly 
seen, one can read off this figure the positions of first energy levels for the 
problem considered. In the literature many different methods of resummation 
of trace formulas were discussed (see e.g. JHl and references therein). 
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Fig. 1. The trace formula for the rectangular billiard with periodic boundary 
conditions calculated by taking into account 250 different periodic orbits. Dotted 
lines indicate the position of exact energy levels. 

2 Billiards on Constant Negative Curvature Surfaces 

The crucial point in the second method of derivation of the trace formula for 
the rectangular billiard with periodic boundary conditions was a representa- 
tion of the exact Green function as a sum of a free Green function over all 
images of the initial point. This method of images can be applied for any 
problem which corresponds to a factorization of a space over the action of a 
discrete group. In the Euclidean plane (i.e. the space of zero curvature) there 
exist only a few discrete groups. Much more different discrete groups are pos- 
sible in the constant negative curvature (hyperbolic) space. Correspondingly, 
one can derive the trace formula (called the Selberg trace formula) for all 
hyperbolic surfaces generated by discrete groups. 

The exposition of this Section follows closely HU]. In Sect . 12 . ll hvperbolic 
geometry is non-formally discussed. The important fact is that on hyperbolic 
plane there exist an infinite number of discrete groups (see e.g. 02]) • Their 
properties are mentioned in Sect. 12.21 In Sect. 12.31 the classical mechanics on 
hyperbolic surfaces is considered and in Sect. l2.4l the notion of quantum prob- 
lems on such surfaces is introduced. The construction of the Selberg trace 
formula for hyperbolic surfaces generated by discrete groups consists of two 
steps. The first is the explicit calculation of the free hyperbolic Green func- 
tion performed in Sect. 12.51 The second step includes the summation over all 
group transformations. In Sect. 12.61 it is demonstrated that the identity group 
element gives the mean density of states. Other group elements contribute to 
the oscillating part of the level density and correspond to classical periodic 
orbits for the motion on systems considered. The relation between group ele- 
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ments and periodic orbits is not unique. All conjugated matrices correspond 
to one periodic orbit. The summation over classes of conjugated elements is 
done in Sect. 12.71 Performing necessary integrations in Sect. 12.81 one gets the 
famous Selberg trace formula. Using this formula in Sect. 12.91 we compute the 
asymptotic density of periodic orbits for discrete groups. In Sect. 12.101 the 
construction of the Selberg zeta function is presented. The importance of this 
function follows from the fact that its non-trivial zeros coincide with eigenval- 
ues of the Laplace-Beltrami operator automorphic with respect to a discrete 
group (see Sect. 12.1111 . Though the Selberg zeta function is defined formally 
only in a part of the complex plan, it obeys a functional equation fScct. in^ 
which permits the analytical continuation to the whole complex plane. 

2.1 Hyperbolic Geometry 

The standard representation of the constant negative curvature space is the 
Poincare upper half plane model {x, y) with y > (see e.g. [7] and 021) with 
the following metric form 

ds^ = l(dx2+d2/2) . 

yZ 

The geodesic in this space (= the straight line) connecting two points is the 
arc of circle perpendicular to the abscissa axis which passes through these 
points (see Fig. ISJ. The distance d(x, y) between two points x = [xi.yi) and 



y 




Fig. 2. The Poincare model of constant negative curvature space. Solid line indicates 
the geodesic passing through points A and B. 



y = (a^2, 2/2) is defined as the length of the geodesic connecting these points. 
Explicitly 
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coshd(x,y) = 1 + i-^--^r + iy^-y^r = 1 + (22) 

22/12/2 2Im zilm Z2 

where in the last equation one combined coordinates (x, y) into a complex 
number z = x + iy. 

In the Euclidean plane the distance between two points remains invariant 
under 3-parameter group of rotations and translations. For constant negative 
curvature space the distance (|22(l is invariant under fractional transformations 

with real parameters a, b, c, d. This invariance follows from the following rela- 
tions 

, , azi +b az2 + b , , , ^ Zi - Z2 
Z\- Z2 = — — = [ad - be) 



and 



ezi + d CZ2 + d {ezi + d) {cz2 + d) 



y' ^ -{z' ~ z'*) = (ad-be)- ^ 



2z ' ' '\ez + d\^ ■ 

Substituting these expressions to (|22|l one concludes that the distance between 
two transformed points z[, z'^ is the same as between initial points zi, 22- 

As fractional transformations are not changed under the multiplication of 
all elements a, 6, c, d by a real factor, one can normalize them by the require- 
ment 

ad—be=\. 

In this case the distance preserving transformations are described by 2 x 2 
matrices with real elements and unit determinant 

^, I , and det q = ad — be ^ 1 . 
e d J ' 

It is easy to check that the result of two successive fractional transformations 
H23(l corresponds to the usual multiplication of the corresponding matrices. 

The collection of all such matrices forms a group called the projective spe- 
cial linear group over reals and it is denoted by PSL(2,IR). 'Linear' in the name 
means that it is a matrix group, 'special' indicates that the determinant equals 
1, and 'projective' here has to remind that fractional transformations (|23l) are 
not changed when all elements are multiplied by ±1 which is equivalent that 
two matrices ±1 corresponds to the identity element of the group. 

The free classical motion on the constant negative curvature surface is 
defined as the motion along geodesies (i.e. circles perpendicular to the abscissa 
axis). The measure invariant under fractional transformations is the following 
differential form 

dM=^. (24) 

y 
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This measure is invariant in the sense that if two regions, D and D' , are 
related by a transformation H23|) . D' = g{D), measures of these two regions 
are equal, n{D') = fJ,{D). 

The operator invariant with respect to distance preserving transformations 
H23() is called the Laplace-Beltrami operator and it has the following form 



for any fractional transformation g{z). 

Practically all notions used for the Euclidean space can be translated to 
the constant negative curvature case (see e.g. [T]). 

2.2 Discrete groups 

A rectangle (a torus) considered in Sect. ^ was the result of the factorization 
of the free motion on the plane by a discrete group of translations (|20|l . Ex- 
actly in the same way one can construct a finite constant negative surface by 
factorizing the upper half plane by the action of a discrete group € PSL{2,]R). 

A group is discrete if (roughly speaking) there is a finite vicinity of every 
point of our space such that the results of all the group transformations (except 
the identity) lie outside this vicinity. The images of a point cannot approach 
each other too close. 



The group of transformation of the unit circle into itself. The group consists 
of all transformations of the following type 



where a is a constant and n is an integer. If a is a rational number a — M/N , 
g{n) can take only a finite number of values (s'(n))^ = 1 and the corresponding 
group is discrete. But if a is an irrational number, the images of any point 
cover the whole circle uniformly and the group is not discrete. 



Modular Group 

Mathematical fact: in the upper half plane there exists an infinite number of 
discrete groups (see e.g. 021 )■ As an example let us consider the group of 2 x 2 
integer matrices with unit determinant 




(25) 



Its invariance means that 



Als/(.9(^))-5(Alb/(^)) 



Example 



z — > g{n)z, and g{n) = exp(27riQ;n) , 
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TO, n, fc, I are integers and ml — nk = 1 . 

This is evidently a group. It is called the modular group PSL(2,K) {TZi means 
integers) and it is one of the most investigated groups in mathematics. 

This group is generated by the translation T : z ^ z-\-l and the inversion 
S : z ^ —1/ z (see e.g. 02) which are represented by the following matrices 

^- (oi) ' ^- (-lo) • 

These matrices obey defining relations 

5*2 = -1 , [STf = 1 

and are generators in the sense that any modular group matrix can be repre- 
sented as a product of a certain sequence of matrices corresponding to S and 
T. 

Fundamental Region 

Similarly to the statement that the rectangular billiard is a fundamental do- 
main of integer translations, one can construct a fundamental domain for any 
discrete group. 

By definition the fundamental domain of a group is defined as a region on 
the upper half plane such that (i) for all points outside the fundamental do- 
main there exists a group transformation that puts it to fundamental domain 
and (ii) no two points inside the fundamental domain are connected by group 
transformations . 

The fundamental domain for the modular group is presented in Fig. |3| In 
general, the fundamental region of a discrete group has a shape of a polygon 
built from segments of geodesies. Group generators identify corresponding 
sides of the polygon. 




2.3 Classical Mechanics 

Assume that we have a discrete group G with corresponding matrices M S 
G e PSL(2,]R) 

The factorization over action of the group means that points z and z' where 

, az + b 



cz + d 



(26) 



are identified i.e. they are considered as one point. The classical motion on the 
resulting surface is the motion (with unit velocity) on geodesies (semi-circles 
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-1-1/2 1/2 1 



Fig. 3. Fundamental domain of the modular group. The indicated parts are identi- 
fied by the corresponding generators 



perpendicular to the real axis) inside the fundamental domain but when a 
trajectory hits a boundary it reappears from the opposite side as prescribed 
by boundary identifications. 

For each hyperbolic matrix M e G with |Tr Af | > 2 one can associate 
a periodic orbit defined as a geodesies which remains invariant under the 
corresponding transformation. The equation of such invariant geodesies has 
the form 



c{x'^ + y^) + {d - a)x - b = 



(27) 



This equation is the only function which has the following property. If z = 
X + iy belongs to this curve then 



cz + d 



also belongs to it. 

The length of the periodic orbit is the distance along these geodesies be- 
tween a point and its image. Let z' as above be the result of transformation 
l|26(l then the distance between z and z' is 



coshlp — 1 



But y' = y/\cz + dp and 

az + b c{x + iy)'^ — (d — a){x + iy) — b 
cz + d cz + d 



—2cy + i{d — a + 2cx) 



Here we have used the fact that point z belongs to the periodic orbit (i.e. its 
coordinates obey H27|l 'l. Therefore 
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cosh Ip = I+2I ^ '^cy + i{d — a + 2acx) \ ^ = 

= l + i[46c+(d-a)2] = i(a + rf)2_i . 

Notice that the length of periodic orbit does not depend on an initial point 
and is a function only of the trace of the corresponding matrix. 
Finally one gets 

2cosh^ = |Tr M| . (28) 

Periodic orbits are defined only for hyperbolic matrices with |Tr M\ > 2. For 
discrete groups only a finite number of elliptic matrices with |Tr Af | < 2 can 
exist (see 

To each hyperbolic group matrix one can associate only one periodic orbit 
but each periodic orbit corresponds to infinite many group matrices. It is due 
to the fact that z and g{z) for any group transformation has to be considered 
as one point. Therefore all matrices in the form 

SMS-^ 

for all e G gives one periodic orbit. These matrices form a class of conju- 
gated matrices and periodic orbits of the classical motion are in one-to-one 
correspondence with classes of conjugated matrices. 

2.4 Quantum Problem 

The natural 'quantum' problem on hyperbolic plane consists in considering 
the same equation as in ^ but with the substitution of the invariant Laplace- 
Beltrami operator H25|) instead of the usual Laplace operator 

for the class of functions invariant (= autoniorphic) with respect to a given 
discrete group G 

where z' = x' + iy' is connected with z — x + iy hy group transformations 

, az + b 

^ = ; • 

cz + a 

It is easy to check that the Laplace-Beltrami operator H25|l is self-adjoint with 
respect to the invariant measure (|24|l . i.e. 

j •f*(A*')d^ = j {M'*)^AiJL 

and all eigenvalues En are real and £'„ > 0. 
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2.5 Construction of the Green Function 

As in the case of plain rectangular billiards the construction of the Green 
function requires two main steps. 

• The computation of the exact Green function for the free motion on the 
whole upper half plane. 

• The summation of the free Green function over all images of the initial 
point under group transformations. 

The free hyperbolic Green function obeys the equation 
(ALB+i?)G(,°\x,x') = <5(x-x') 
and should depend only on the (hyperbolic) distance between points x, x' 



y = cosh(i(x, x') = 1 + 



After simple calculations one gets that G{y) with y obeys the equation 
for the Legendre functions (see e.g. Vol.1, Sect. 3) 



— ~ 2y— 
dy^ dy 



(l-2/')33-~2y^ + K^ + l)G = 

where 



E=\ + k^ = -l{l + l) 
4 

and 

l^---ik . 
2 

As for the plane case the required solution of the above equation should grow 
as e*'""'* when d ^ oo and should behave like ln(i/27r when d ^ Q. From [62\ . 
Vol.1, Sect. 3 it follows that 

G'e'(x,x') = -^Q_i_ij;(coshd(x,x')) . 

Here (5_i_jj.(cosh(i) is the Legendre function of the second kind with the 
integral representation [S^], Vol. 1 (3.7.4) 



1 f° 

Q_i_ifc(coshd) = -yiy^ 

and the following asymptotics 



V cosh r — cosh d 



_ jj, (cosh d) '^—^ — log d 



and 
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The automorphic Green function is the sum over aU images of one of the 
points 

G£;(x,x')=^4°^(x,<7(x')) 
9 

where the summation is performed over all group transformations. 

2.6 Density of State 

Using the standard formula H17(l 

d{E) = -- [ Im Gij(x,x)d/i 

71" Jd 

one gets the expression for the density of states as the sum over all group 
elements 

dxdy I f°° sin krdr \ 



a 

Mean Density of States 



y \Jdiz,giz)) i/coshr - cosh d{z,g{z)) 



The mean density of states corresponds to the identity element of our group. 
In this case g{z) — z and d{z,g{z)) = 0. Therefore 

diE) = ^f^r -^^dr 
2y27r2 Jd y^ Ja Vcoshr- 1 

li(D) f°° sinfcr , 
' dr 



dxdy 

D y'^ 

is the (hyperbolic) area of the fundamental domain. 
The last integral is 

sinfcr 

d?' ~ TT tanh nk 



(27r)2 sinh(r/2) 

where 



/g sinh(r/2) 
and the mean density of states takes the form 

d(E) = ^^^tanhTrfc . 
47r 

When fc — > cx) it tends to ^{D)/Att as for the plane case. 
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2.7 Conjugated Classes 

The most tedious step is the computation of the contribution from non-trivial 
fractional transformations. 

Let us divide all group matrices into classes of conjugated elements. It 
means that all matrices having the form 

g' = SgS-^ 

where S belong to the group are considered as forming one class. 

Two classes either have no common elements or coincide. This statement 
is a consequence of the fact that if 

SigiSi ^ = S2g2S2 ^ 

then g2 = SsgiS^^ where S3 = S'f ^5*2. Therefore (?2 belongs to the same class 
as gi and group matrices are splitted into classes of mutually non-conjugated 
elements. 

The summation over group elements can be rewritten as the double sum 
over classes of conjugated elements and the elements in each class. Let 5 be a 
representative of a class. Then the summation over elements in this class is 

J2 I f[z,SgS-\z))di, 
s 

and the summation is performed over all group matrices S provided there is 
no double counting in the sum. The latter means that matrices S should be 
such that they do not contain matrices for which 

SigS^^ = S2gS2^ 

or the matrix S3 = Si^S2 commutes with matrix g 

Ssg = gSs ■ 

Denote the set of matrices commiiting with g by Sg . They form a subgroup 
of the initial group G as their products also commute with g. To ensure 
the unique decomposition of group matrices into non-overlapping classes of 
conjugated elements the summation should be performed over matrices S such 
that no two of them can be represented as 

S2 = sSi 

and s belongs to Sg. This is equivalent to the statement that wo sum over all 
matrices but the matrices sS are considered as one matrix. It means that we 
factorize the group over Sg and consider the group G/Sg. 

As the distance is invariant under simultaneous transformations of both 
coordinates 
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d{z,z') = d{Siz),Siz')) 

one has 

diz,g{z)) = d{S{z),Sg{z)) = diy,SgS-\y)) 

where y = S{z). 

These relations give 

/ f{d{y,SgS-\y)))di^= [ f{z,g{z))A,i 
Jd Js-''(d) 

and the last integral is taken over the image of the fundamental domain under 
the transformation S~^. Therefore 

f{d{y,SgS-\y)))df, = Y, [ mz,g{z)))dti . 

For different S images S~^{D) are different and do not overlap. The integrand 
does not depend on S and 

^/ f{d{y,SgS-\y)))di,= I f{d{z,g{z)))d^i 
s Jd JDg 

where 

s 

The sum of all images S~^{D) will cover the whole upper half plane but we 

have to sum not over all S but only over S factorizcd by the action the group 
of matrices commuting with a fixed matrix g. Therefore the sum will be a 
smaller region. 

Any matrix g can be written as a power of a primitive element 

5 = 5o" 

and it is (almost) evident that matrices commuting with g are precisely the 
group of matrices generated by go- This is a cyclic abelian group consisting of 
all (positive, negative, and zero) powers of 170 

Sg=g^, m = 0,±l,±2,... 

and as a discrete group it has a fundamental domain FDg. 
Therefore 

J2 f fidiy,SgS-'iy)))d^i= [ f{d{z,g{z)))dtx. 
seG/s, JD 

In the left hand side the integration is taken over the fundamental domain of 
the whole group G and the summation is done over matrices from G factorized 
by the subgroup Sg of matrices which commutes with a fixed matrix g. In the 
right hand side there is no summation but the integration is performed over 
the (large) fundamental domain of the subgroup Sg. 
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2.8 Selberg Trace Formula 

We have demonstrated that the density of states of the hyperbohc Laplace- 
Beltrami operator automorphic over a discrete group can be represented as 

d{E) =d{E)+J2dgiE) 
a 

where 

d^^E) = [ d^ r — 

2\/2n'^ JpDg Jd(z,g(z)) vcoshr - coshd(z,5(z)) 

and the summation is performed over classes of conjugated matrices. 

Let us consider the case of hyperbohc matrices 5 = 5™ (i-^- |Tr .go| > 2). 
By a suitable matrix B such matrix can be transform to the diagonal form 



Ao 
A„-i 



For hyperbolic matrices Aq is real and |Ao| > 1. By the same transformation 
the matrix g will be transformed to 



A 

A-i 



BgB-^ 
and A = AJJ'. 

Assume that g is in the diagonal form. Then g{z) = A^z and 

coshd(z,5(z)) = 1 + ^ ^^^^^ ^ ' . 

Because Ao is real the transformation z' — AqZ gives y' = AqJ/ and the funda- 
mental domain of Sg = Aq™z has the form of a horizontal strip 1 < y < Aq 
indicated in Fig. 0] Now 



oo /"A 



Introducing a new variable ^ — xy one gets 



y J~oo \ ^ 
A^ 



After the substitution 
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Fig. 4. Fundamental domain of multiplication group 



A2 

one obtains 

lnA2 f°° F{u) 



where 



dg{E) = / "^ du 

-^"0 V w - ""0 



(A2-1)2 1 



A2 ' A2 

The variable u is connected with the distance by coshd = l + u/2 and the 
function F(coshrf) has the form 

i^(coshd) = — [ sin kr ^ 

2\/27r2 Jd vcoshr — coshd 

Introduce a variable r connected with r as m is connected with d 

r dr 1 



cosh r = 1 



2 ' dr + 4r 



It gives 



and 



where 



p , , 1 p sin fer(T) 
-F(u) = -— r / dr 
27r2y„ ^(r-u)(r2 + 4r) 



Vu^^ Ju ^J{t-u){t'^ +4r) 
Changing the order of integration one obtains 
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sin kr{T) dt 



Jw Vt^+^T Ju 



The last integral is a half of the residue at infinity 

du 



/ 



= TT 



\/{u — w)(t — u) 
and 

/" Sin JcT" ( T) f 71" 

fiw) = TT =dT = TT sin(fcr)dr = — cos fcL . 

' L V^^T4^ J I, k " 

Here Ip is the minimal value of r corresponding to uq 

coshZ, = l + f = l + i(A^ + 1-2) = i(A+ 1)^-1 

or 

2coshZp = A + - = Tr 5 
A 

i.e. Ip is the length of periodic orbit associated with the matrix g. 
Therefore 

InA^ ^(0) 

dniE) = ° cos klr, = — ; — ^ , , , cos kL 

' 27rfcVA + A-i-2 47r/csinh/p/2 " 

where l^p^ is the length of the primitive periodic orbit associated with go- 

Combining all terms together one finds that the eigenvalues density of the 
Laplace-Beltrami operator automorphic with respect to a discrete group with 
only hyperbolic matrices has the form 

m(^) , V- ^ cos{knlp) 



d{E) = ^^^tanh7rfc+ ^ 



Att ^ ATvk sinh(nL/2) 

p.p.o. n=l \ fi / 

The oscillating part of the density is given by the double sum. The first sum- 
mation is done over all primitive periodic orbits (p.p.o.) and the second sum is 
performed over all repetitions of these orbits. Here k is the momentum related 
with the energy by = fc^ + 1/4. 

To obtain mathematically sound formula and to avoid problems with con- 
vergence it is common to multiply both parts of the above equality by a test 
function h{k) and to integrate over dE = 2kdk. To assume the convergence 
the test function h{r) should have the following properties 

• The function h{r) is a function analytical in the region |Im r\ < 1/2 + 6 
with certain 5 > 0. 

• h{—r) = h{r). 

• \h{r)\<A{l + \r\)-^-^. 
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The left hand side of the above equation is 

J d{E)h{k)AE ^^5{E~ E„)h{k)dE = ^ /i(fc„) . 

n n 

In the right hand side one obtains 

/ h{k)^^kdk r h{k)e-"''dk . 

J 27rfc 27r 

The final formula takes the form 



V/i(fc„) = ^^^/ fc/i(fc) tanh(7rfc)dfc 
„ 27r J„oc 



OO 
OO 



+ j:ipT.j-^f(-i,) (29) 

p.p.o. n=l \ fi / 

where fc„ is related with eigenvalue E„ as follows 

En = kl + - 
and g{l) is the Fourier transform of h{k) 

gil) = / /j(fc)e-'«dfc . 

This is the famous Selberg trace formula. It connects eigenvalues of the 
Laplace-Beltrami operator for functions automorphic with respect to a dis- 
crete group having only hyperbolic elements with classical periodic orbits. 

2.9 Density of Periodic Orbits 

To find the density of periodic orbits for a discrete group let us choose the 
test function h{r) in as 

h{r) = e-^'-'+i/^)^ ^ e-^^ 
with a parameter T > 0. fts Fourier transforms is 

g{u) = — / h{k)e-"'^dk = ^— e-"'/4^ . 

In the left hand side of the Selberg trace formula one obtains 



n E„>0 
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where we take into account that for any discrete group there is one zero 
eigenvalue corresponding to a constant eigenfunction. Therefore when T ^ oo 
the above sum tends to one 

^e-^-^^-^1. 

n 

One can easily check that in the right hand side of H29|) the contribution of 
the smooth part of the density goes to zero at large T and the contribution 
of periodic orbits is important only for primitive periodic orbits with n = I. 
The latter is 

T/4 a—T/4 foo 



2%AfT ^ " 2\/7rT 

where p{l) is the density of periodic orbits. Hence the Selberg trace formula 
states that 

-T/4 



hm / le-'-/'^'''^p{l)dl = 1 



Assume that p{l) — he°'^ /I with certain constants a and b. Then from the 
above limit it follows that a = 6 = 1 which demonstrates that the density of 
periodic orbits for a discrete group increases exponentially with the length 

e' 

P{1) = T ■ 



2.10 Selberg Zeta Function 

Amongst many applications of the Selberg trace formula let us consider the 
construction of the Selberg zeta function. 
Choose as test function /i(fc) the function 

Its Fourier transform is 

' 2a 2(3 
The Selberg trace formula gives 

1 



^ ( i.2 I 



p{D) 



k tanh Trfc — — — — dk 



2-K \k^ + a^ k^+(3^' 

Ir 



^ ^ 2sinhnL/2 V 2a 23 

p.p.o. n—1 f/\ 
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The Selberg zeta function is defined as the following formal product 



(30) 



p.p.o m=0 



One has 

1 dZ 
Z ds 



Z—/ \ _ Q—lp{s+m) 
p.p.o. m— p.p.o. n=l m=0 



>.p.o. ?i=l m=0 

, f^—lpUS ^ 1 

V Z V — = V Z V „-i„n(s-l/2) 

4-'„ ^^l-e-'p" ^ ^^2sinhnL/2 



p.p.o. n=l " p.p.o. n=l 

Choose a = s — 1/2 and (3 = s' — 1/2 then 
1 1 



E 



+ (s - 1/2)2 ^2 + (5/ _ 1/2)2 



47r 



fc tanh Trfc 



1 



fc2 + (s- 1/2)2 A;2 + (s'- 1/2)2 



dA; + 



+- 



1 Z{s)' 1 Z(s' 



2s - 1 Z{s) 2s' - 1 Z(s') 
The integral in the right hand side can be computed by the residues 

1 1 



j A; tanh TT A; ^ 



fc2 + (s- 1/2)2 A;2 + (s- 1/2)2 



dfc = /(s) - /(s') 



where /(s) is the sum over residues from one pole k = i{s — 1/2) and from 
poles fc„ = i(n + 1/2) of tanhTrfc 



f{s) = 27ri 



tanh[i7r(s - 1/2)] + - V 



n + l/2 



00 ^ cx> ^ 

TT cot — h . 



s — n 5 + n 

n=l n=l 



But 



00 -j^ 00 

TT cot TTS = > h > 

^-^ x — n ^-^ -I- 



, s — n — ' s + n 

n=l n=l 



therefore 



00 ^ 
/(«)=2E7^ 



s + n 



Using these relations one gets the identity valid for all values of s and s' 
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1 Z'(s) _ 1 Z'(s') ix{D) g 1 



2s - 1 Z{s) 2s' ~ 1 Z(s') 27r ^f:^ V « + « s' + n 



^ ? (^^ + (-^ - 1/2)2 ^.2 + _ 1/2)2) • (31) 

The right hand side of this identity has poles at s = 1/2 + ifc„ and s — —n. 
The same poles have to be present in the left hand side. If 

Z'(s) i^k 



Z{s) s-Sk 

then 

^i^) ^ ^ SkY'' when s Sk ■ 

When I'k > (resp. z/fc < 0) point Sk is a zero (resp. a pole) of the Selberg 
zeta function Z{s). 



2.11 Zeros of the Selberg Zeta Function 

Combining all poles one concludes that the Selberg zeta function for a group 
with only hyperbolic elements have two different sets of zero. The first consists 
of non-trivial zeros 

s = l/2±ifc„, 

coming from eigenvalues of the Laplace-Beltrami operator for automorphic 
functions. The second set includes a zero from E = eigenvalue and zeros 
from the smooth term. These zeros are called trivial zeros and they are located 
at points 

s = — m (to = 1,2,... 

with multiplicity Vrn = (2to + l)/K(Z3)/27r, at point s = with multiplicity 
vq — fi{D)/2'K and a single zero at s = 1. These multiplicities are integers 
because the area of a compact fundamental domain /i(-D) = 47r(g — 1) where 
g is the genus of the surface. 

The structure of these zeros is presented schematically at Fig. 



2.12 Functional Equation 

The infinite product defining the Selberg zeta function (|30|l converges only 
when Re s > 1/2. Nevertheless the Selberg zeta function can be analytically 
continued to the whole complex plane s with the aid of 1)31(1 . 

Put s' = 1 — s in ((31|l . The sum over eigenvalues cancels and f{s) — /(I — 
s) — 27rcot7rs. Therefore 

1 fZ'is) Z'{l-s)\ fi{D) 
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non-trivial zeros 



-4 -3 -2 -1 



Fig. 5. Zeros of the Selberg zeta function 

which is equivalent to the following relation (called functional equation) 

Z{s) = ip{s)Z{l- s) (32) 

where 



^^-^,{D){s-l-) cot ITS 
ip{s) 2 



and 93(1/2) = 1. 
Explicitly 



(/5(s) = exp n{D) 



-1/2 



utamrudu 



Therefore if one knows the Selberg zeta function when Re s > 1 (|32|) gives its 
continuation to the mirror region Re s < 0. 



3 Trace Formulas for Integrable Dynamical Systems 

A /-dimensional system is called integrable if its classical Hamiltonian can be 
written as a function of action variables only 

H{I)=H{Iu...,If). 

In this representation the classical equations of motion take especially simple 
form 

t OH dH 

The semiclassical quantization consists of fixing the values of the action vari- 
ables 
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I3 = Hn, + f ) 

where Uj are integers and fij are called the Maslov indices. 

In this approximation eigenvalues of energy of the system are a function 
of these integers 

The eigenvalue density is the sum over all integers rij 
d(ii;)=^<5(i?-i/(n(n+iM))- 

n 

Using the Poisson summation formula one transforms this expression as 
follows 

d{E) =J2 [ e'^'^'SiE - H{h{n + i/x))dn = 

e^'"^''/' J e^'^'Ni/;. J _ H{I))dI (33) 
where the summation is taken over / integers Nj. 
3.1 Smooth Part of the Density 

The term with N = in H33|) corresponds to the smooth part of the density 



diE) = ^ J S{E-H{I))dI 



As dld<p is the canonical invariant, dldtp = dpdq where p and q are the 
momenta and coordinates and, because / d(p ~ (27r)-'', the formula for the 
smooth part of the level density can be rewritten in the Thomas-Fermi form 

diE)= J SiE-Hip,<i))-^^. (34) 

The usual interpretation of this formula is that each quantum state occu- 
pies {2Tih)f volume on the constant energy surface. For general systems 134|l 
represents the leading term of the expansion of the smooth part of the level 
density when h ^ Q. Other terms can be found e.g. in (51 . See also for 
the resummation of such series for certain models. 
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3.2 Oscillating Part of the Density 

In the semiclassical approximation fi ^ terms with N 7^ in H33|) can be 
calculated by the saddle point method. Our derivation differs slightly from 
the one given in [^. First it is convenient to represent (5- function as follows 



00 

e 

00 



Then 



N — C30 J 



where the effective action, S{l,a), is 

S'(I, a) ^ 27rNI + a{E - H{I)) . 

The integration over I and a can be performed by the saddle point method. 
The saddle point values, I^p and asp, are determined from equations 

— = £; - H{Isp) = , — = 2^N - aspUJsp = . 

The first equation shows that in the leading approximation I^j, belongs to the 
constant energy surface and the second equation selects special values of I^p 
for which frequencies coj are commensurable 

t^sp = N . 

^sp 

Together the saddle point conditions demonstrate that in the limit ?i — > the 
dominant contribution to the term with fixed integer vector N comes from 
the classical periodic orbit with period 

Tp — Ckg.p 

and the saddle point action coincides with the classical action along this tra- 
jectory 

Ssp = 27rNIsp . 

To compute remaining integrals it is necessary to expand the full action up 
to quadratic terms on deviations from the saddle point values. One has 



T 



S{Isp + SI, as.p + 5a) = S^p + -^(SI^Hi^SIj) - 5a{ujj5I. 



where the summation over repeating indexes is assumed. Hij is the matrix of 
the second derivatives of the Hamiltonian computed at the saddle point 
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The following steps are straightforward 



dSIdSaexp [ ^S{I,a) 



I 

//2 



V^?'/ ^|deti/„|(a;,i/,^i.;,^ V?^ 4 
(27r)(/-i)/2?i(/+i)/ 



2 

sp 



exp ( — + -7Tf3 



where Qij = H~j^ det H called the co-matrix of Hij is the determinant ob- 
tained from Hij by omitting the i-th row and the j-th column. The phase /? 
is the signature of Hij minus the sign of {ujH^^uj). 

The final expression for the oscillating part of the level density of an inte- 
grable system with a Hamiltonian H{T) is 

d(-) (E) =Y^P^ exp (i| - i JNm + i J/?) 

where Sp = 27rNI is the action over a classical periodic orbit with fixed 
winding numbers and 



2^ X (/-3)/2 



The summation over integer vectors N is equivalent to the summation over 
all classical periodic orbit families of the system. 



4 Trace Formula for Chaotic Systems 

To compute the eigenvalue density for a chaotic system one has to start with 
general expression itTTIl 

d{E) = -\ j Im G£;(x,x)dx 

which relates the quantum density with the Green function of the system, 
GB(x,y), obeying the Schroedinger equation with a (5-function term in the 
right hand side 
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For concreteness let us consider the usual case 

H = -h^A + V{yi) . 

The exact Green function can be computed exactly only in very limited cases. 
For generic systems the best which can be achieved is the calculation of the 
Green function in the semiclassical limit h —i- 0. 

4.1 Semiclassical Green Function 

Let us try to obey the Schroedinger equation in the following form (see 

G£(x,y) = ^(x,y)e'^(-^y)/'' (35) 

where the prefactor ^(x.y) can be expanded into a power series of h. 

Separating the real and imaginary parts of the Schroedinger equation one 
gets two equations 

{E - (VS')2 - y(x)) + h^AA = 

and 

2V5VA + ASA = . 

In the leading order in h the first equation reduces to the Hamilton- Jacobi 
equation for the classical action S'(x, y) 

E = (VS')2 + V{x) . 

It is well known that the solution of this equation can be obtained in the 
following way. 

Find the solution of the usual classical equations of motion 

dV 

with energy E which starts at a fixed point y and ends at a point x. Then 



5(x,y) = 




where p is the momentum and the integral is taken over this trajectory. 

Instead of proving this fact we illustrate it on an example of the free 
motion. The free motion equations x = have a general solution 

X = + y 

with a fixed vector k. One has 
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t 



and the conservation of energy |k|' 



E determines the time of motion 

<;-y| 



E 



Therefore 



5(x,y) = Vl|x-y| 



which, evidently, is the solution of the free Hamilton-Jacobi equation. 
The next order equation 

2V5'Vyl + ASA ^ 

is equivalent to the conservation of current. Indeed, for the semiclassical wave 
function lP5]l 



2V 



and 



VJ = ^(2VylVS' + AAS) = . 
The solution of the above transport equation has the form 



^(x,y) 



(27r?i)(/+i)/2 



1 



■det 



dtj I dt 



i±Uif± 



1/2 



where ti± and tf± are coordinates perpendicular to the trajectory in the 
initial, y, and final, x, points respectively and ki, kf are the initial and final 
momenta. The derivation of this formula can be found e.g. in . The overall 
prefactor in this formula can be fixed by comparing with the asymptotics of 
the free Green function ((TUl at large distances. 

The final formula for the semiclassical Green function takes the form 



G£(x,y) 



E 



classical 
trajectories 



(27ra)(/+i)/2 



1 



hi kj 



■ dot - 



dti^dt 



1/2 



X exp ( ^5c/(x,y) - ^tt^ 



where the sum is taken over all classical trajectories with energy E which 
connect points y and x. /j, is the Maslov index which, roughly speaking, counts 
the number of points along the trajectory where semiclassical approximation 
cannot be applied. 
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4.2 Gutzwiller Trace Formula 

The knowledge of the Green function permits the calculation of the density 
of eigenstates by the usual formula (fT7|) 

d{E) = -^ J Im Ge (x, x)dx . 

The Green function G£;(x,y) at points x and y very close to each other 
has two different contributions (see Fig. EJ. The first comes from very short 



Small action trajectory 




Large actiontrajectory 



Fig. 6. Small and large action contributions to the Green function for nearby points 



trajectories where semiclassical approximation cannot, in general, be applied. 
The second is related with long trajectories. The first contribution can be 
computed by using the Thomas-Fermi (local) approximation for the Green 
function. In this approximation one uses the local formula (cf. (llGfl l 



GEi^.yV-^ j 



dp c*p(x-y)/?i 
{2Tiny {E- H{p,x) 4 



Therefore 



Im Gs(x,x) = -TT J j^^S{E - H{p,x 



)) 



and the smooth part of the level density in the leading approximation equals 
the phase-space volume of the constant energy surface divided by {2'Khy 



The contributions from long classical trajectories with finite actions corre- 
sponds to the oscillating part of the density and can be calculated using the 
semiclassical approximation of the Green function (|35|l . 
One has 
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(£;) = -- Im y"A(x,x)e'^(^''')/''dx 



classical 
trajectories 



When h the integration can be performed in the saddle point approxi- 
mation. The saddles are solutions of the equation 



dS{^,y) , dS{x,y) 



dx dy 



= 



But 



a^(x,y) , 95(x,y) 



dx ' ^ ' dy 

where kj and are the momenta in the final and initial points respectively. 

Hence the saddle point equations select special classical orbits which start 
and end in the same point with the same momentum. It means that the saddles 
are classical periodic orbits of the system and 

Sgp — Sp . 

To calculate the integral around one particular periodic orbit it is convenient 
to split the integration over the whole space to one integration along the orbit 
and (/ — I) integrations in directions perpendicular to the orbit. For simplicity 
we consider the two-dimensional case. 

The change of the action when a point is at the distance y from the periodic 
orbit is 



1 ^ d^S{y,y) ^ 
°^ = — 7^ — la=o 



(9y2 

where S{y, y) is the classical action for a classical orbit in a vicinity of the 
periodic orbit (see Fig. [TJ. To compute such derivatives it is useful to use 




claisTcal orbit 

Fig. 7. A periodic orbit and a closed classical orbit in its vicinity 

the monodromy matrix, my, which relates initial and final coordinates and 
momenta in a vicinity of periodic orbit in the linear approximation 
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Spf J \ m2i m22 J \ Spi 

As the classical motion preserves the canonical invariant dpdq it follows that 
detM= 1. 



One has 



But 



Therefore 



Syf = muSyi + mi26pi , 

5pf = TO21% + m22Spi . 

dS dS 
ovi oyf 



From comparison of these two expression one obtains the expressions of the 
second derivatives of the action through monodromy matrix elements 

d'^S 1_ d^_rnn d^S _ m22 

dyidyf mi2 ' dy"? mi2 ' dyj mu ' 

Substituting these expressions to the contribTition to the trace formula from 
one periodic orbit one gets (in two dimensions) 

1 /■ I i_i/2 /if, , .m,ii+m22-2 on, da; 



k{x) 



where x and y are respectively coordinates parallel and perpendicular to the 
trajectory. 

Computing the resulting integrals one obtains 

^ 7r7i i/lmii +m22 - 2| 

where Tp = J dx/k{x) is the geometrical period of the trajectory. 

Finally the Gutzwiller trace formula takes the form (valid in arbitrary 
dimensions) 



T °° 1 

primitive 
periodic 
orbits 



In the derivation of this formula we assumed that all periodic orbits are un- 
stable and Mp is the monodromy matrix for a primitive periodic orbit. 
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5 Riemann Zeta Function 



The trace-like formulas exist not only for dynamical systems but also for the 
Riemann zeta function (and others number-theoretical zeta functions as well) . 

The Riemann zeta function is a function of complex variable s defined as 
follows 

^i')-f:^ = ]li^-p''r' (36) 

n—1 p 

where the product is taken over prime numbers. The second equality (called 
the Euler product) is a consequence of the unique factorization of integers 
into a product of prime numbers. 

This function converges only when Res > 1 but can analytically be con- 
tinued in the whole complex s-plane. 



5.1 Functional Equation 

The possibility of this continuation is connected with the fact that the Rie- 
mann zeta function satisfies the important functional equation 

as) = ^{s)ai - s) (37) 

where 

^(s) = 2%«-isin(^)r(l-s). (38) 

We present one of numerous method of proving this relation (see e.g. |55p. 
When Re s > one has the equality 



oo 



where Fix) is the Gamma function (see e.g. jHU, Vol. 1, Sect. 1). Therefore 
if Re s > 1 

^(^/2)C(g) r sl2^1,j. ^ , 

where iZ^(x) is given by the foUowing series 

oo 
n=l 

Using the Poisson summation formula lOl one obtains 

oo oo 
n— — oo 

which leads to the identity 



/X 

n— — oo 
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2!f(.T) + 1 = ^ (2if(-) + 1 



Hence 



1 

2 Jo Ji 

(^x-'/^-y^+x'/^-^^^{x)dx . 



s — 1 s Jq X 
1 



sis -I] 



The last integral is convergent for all values of s and gives the analytical 
continuation of the Riemann zeta function to the whole complex s-plane, the 
only singularity being the pole at s = 1 with unit residue 

/-I \ ^ 
C(s) - 



s- 1 



(The pole at s = is canceled by the pole of r{s/2) giving ({0) — —1/2.) 

One of important consequences of the above formula of analytical contin- 
uation is that it does not change under the substitution s 1 — s. Therefore 
for all values of s 

or 

Cis) = ipis)C{l - s) 

where 

r[s/2) 

By standard formulas (see e.g. 132. , Vol. 1, 1.2.5, 1.2.15) 

r(x)r(i ~x)^ , r{2x) = 2^^-^Tr-^/^r{x)r{x + 1) 

sm TTX 2 

the last expression can be transformed to (|38|l which proves the functional 
equation (|37|l . 

From the functional equation (|37|) is follows that C(s) has 'trivial' zeros 
at negative even integers (except zero) s — —2, —4, . . . which appear from 
sin(7rs/2) in (jSHJ- All other non-trivial zeros, Ci^n) — 0, are situated in the so- 
called critical strip < Re s < 1. If one denotes these zeros as s„ — 1/2 + i-fn 
then functional equation together with the fact that C(s)* — C(s*) state that 
in general there exit 4 sets of zeros: 7„ , — 7„ , 7* , —7*. 

According to the famous Riemann conjecture (see e.g. |2S]) all nontrivial 
zeros of C(s) lie at the symmetry line Re s = 1/2 or 7„ are all real quantities. 
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Numerical calculations confirms this conjecture for exceptionally large number 
of zeros (see e.g. |47| and the web site of Odlyzko but a mathematical 
proof is still absent. 



5.2 Trace Formula for the Riemann Zeros 

Let us fix a test function h{r) exactly as it was done for the Selberg trace 
formula in Sect. I2.8l i.e. 

• h{r) is a function analytical in the region |Im r\ < 1/2 + 5, 

• h{—r) = h{r), 

• \h{r)\ < A{1 + \r\)-^-^ . 



Denote as in that Section 



1 



+ 00 



= ^ / /i(r)e-'™dr 



OO 



and define 



r' + OO 



Now let us compute the integral 



where the contour of integration is taken over the rectangle — r/< Res<H-?7 
and —T < Im s < T with < rj < S and T — > +oo. Inside this rectangle 
there are poles of C{s)/Cis) coming from non-trivial zeros of the Riemann 
zeta function, s„ = 1/2 + i7„, and the one from the pole of C(s) at s = 1. The 
total contribution from these poles is 

n 

One can check that the limit T oo exists and, consequently, one has the 
identity 



E 



2 27n7i+,,_i^ C(s) 27nJ_,^_i^ ({s) 



Let us substitute in the second integral the functional equation (|37|l with (p{s) 
from H39|l . One has 



as) c(i - s) 2 
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Now all integrals converge and one can move the integration contour till s 
1/2 + i?' with real r. In this manner one obtains 



dsH{s) 

7]—ioo 



T' 




T 


(1) 



r' fi-s 



r \ 2 



1, 1 /■+°° . , - r /I 1 



47ri 

The first term in the right hand side of this equality is due to the appearance 
of the pole of r{s/2) at s = when the integration contour shifted till s = 
1/2 + ir. Also we have used that h{—r) — h{r). 

For terms with the Riemann zeta function one can use the expansion which 
follows from 

' p n=l 

Shifting the integration contour as above (i.e. till s = 1/2 + ir), using that 
g{—u) = g(u), and combining all terms together one gets the following Weil 
explicit formula for the Riemann zeros 

1 /"OO p/ / 1 * \ 

^ '.h,.)^^/^/.W-p(j + ^r)d,. + /,(i) + M-i 

non-trivial 



oo 



5(0) In^- 2 J2 ^^pY.: 

— — g(nlnp) . 



n=l 



Here 7„ are related with non-trivial zeros of the Riemann zeta function, s„, 
as follows 

1 

= - + l-fn ■ 

This formula is an analog of usual trace formulas as it relates zeros of the 
Riemann zeta function defined in a quite complicated manner with prime 
numbers which are a common notion. 

The similarity with dynamical trace formulas is more striking if one as- 
sumes the validity of the Riemann conjecture which states that 7„ are real 
quantities (which in a certain sense can be considered as energy levels of a 
quantum system). In 'semiclassical' limit r — s- oo using the Stirling formula 
(see e.g. Vol. 1, 1.9.4) 

In r(z) '^-tr (z - 1) In z - z + i In 27r 
^ ' 2' 2 

one obtains that the density of Riemann zeros 

d(i?)=^<5(ii;-7n) 
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can be expressed by the following 'physical' trace formula valid at large E 

d{E) = d{E) + d^'^'^HE) 

where 

- \ E 

d{E) = — In h corrections , 

and 

d(-)(i?) = -i5:£^cos(i?nlnp) 

p n—l 

where the summation is performed over all prime numbers. 



5.3 Chaotic Systems and the Riemann Zeta Function 

By comparing the above equations with the trace formulas of chaotic systems 
one observes (see e.g. QSl) a remarkable correspondence between 

different quantities in these trace formulas 

• periodic orbits of chaotic systems <-!■ primes, 

• periodic orbit period Tp ^ \np, 

• convergence properties of both formulas are also quite similar. 

The number of periodic orbits with period less than T for chaotic systems is 
asymptotically 

hT 

NiTp < T) = 

where a constant h is called the topological entropy. 

The number of prime numbers less than x is given by the prime number 
theorem (see e.g. 




As Inp = Tp this expression has the form similar to number of periodic orbits 
of chaotic systems with h = l 

N{Tp < T) = ^ . 

Due to these similarities number-theoretical zeta functions play the role of a 
simple (but by far non-trivial) model of quantum chaos. 

Notice that the overall signs of the oscillating part of trace formulas for 
the Riemann zeta function and dynamical systems are different. According to 
Connes I^OI it may be interpreted as Riemann zeros belong not to a spectrum 
of a certain self-adjoint operator but to an 'absorption' spectrum. Roughly 
speaking it means the following. Let us assume that the spectrum of a 'Rie- 
mann Hamiltonian' is continuous and it covers the whole axis. But exactly 
when eigenvalues equal Riemann zeros corresponding eigenfunctions of this 
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Hamiltonian vanish. Therefore these eigenvalues do not belong to the spec- 
trum and Riemann zeros correspond to such missing points similarly to black 
lines (forming absorption spectra) which are visible when light passes through 
an absorption media. In Connes' approach the 'Riemann Hamiltonian' may be 
very simple (see also ^Sl) but the functional space where it has to be defined 
is extremely intricate. 

6 Summary 

Trace formulas can be constructed for all 'reasonable' systems. They express 
the quantum density of states (and other quantity as well) as a sum over classi- 
cal periodic orbits. All quantities which enter trace formulas can be computed 
within pure classical mechanics. 

Trace formulas consist of two terms 

d{E) = d{E)+d^'"'\E). 

The smooth part of the density, d{E), for all dynamical systems is given by 
the Thomas-Fermi formula (plus corrections if necessary) 

^"(^)-/(W^(^-^(p'-»- 

For integrable systems the oscillating part of the density, d'^°^'^\E), is 

_ / o„ \ (/-3)/2 -, 

^\f^Tj n\mQ,,N,)\ 

where Sp = 27rNI is the action over a classical periodic orbit with fixed 
winding numbers N and Qij is the co-matrix of the matrix of the second 
derivatives of the Hamiltonian. 

For chaotic systems d^°'^'^\E) is represented as a sum over all classical 
periodic orbits 

d(°-)(i<;) IlLy ^ cos (n^ ~ n-fi^] 

^ ' ^ TT^i ^ det Af" - 1 1/2 I n 2^"} 

p.p.o. n— 1 I V p / 1 ^ / 

where Sp is the classical action along a primitive periodic trajectory and Mp 
is its monodromy matrix. 

Usually trace formulas represent the dominant contribution when ?i ^ 0. 
They are exact only in very special cases as for constant negative curvature 
surfaces generated by discrete groups where they coincide with the Selberg 
trace formula. For a group with only hyperbolic elements 



4 



d(E) = ^IL^tanhTTfc 
An 
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where /x(-D) is the area of the fundamental domain of the group and 

Aosc)(T7\ ^ cos{knlp) 
^ ' ^ Amk^ sinh(n?„/2) 

p.p.o. n— 1 \ fi J 

where 1^ are lengths of periodic orbits. 

The formulas similar to trace formulas exist also for number-theoretical 
zeta functions (assuming the generalized Riemann conjecture). In particular, 
for the Riemann zeta function 

and 

prime n— 1 

The principal difficulty of all trace formulas is the divergence of the sums over 
periodic orbits. To obtain a mathematically meaningful formula one considers 
instead of the singular density of states its smoothed version defined as a 
sum over all eigenvalues of a suitable chosen smooth test-function. When 
its Fourier harmonics decrease quickly the resulting formula represent a well 
defined object. 

Suggestions for Further Readings 

• A very detailed account of trace formulas derived by multiple scattering 
method can be found in a series of papers by Balian and Bloch 8 . 

• A short but concise mathematical review of hyperbolic geometry is given 
in 17. 

• Explicit forms of the Selberg trace formula for general discrete groups with 
elliptic and parabolic elements are presented in two volumes of Hejhal's 
monumental work 1391 which contains practically all known information 
about the Selberg trace formula. 

• In PH] one can find a mathematical discussion about different relations 
between number-theoretical zeta functions and dynamical systems. 



2 



Statistical Distribution of Quantum 
Eigenvalues 



Wigner and Dyson in the fifties had proposed to describe comphcated (and 
mostly unknown) Hamiltonian of heavy nuclei by a member of an ensemble 
of random matrices and they argued that the type of this ensemble depends 
only on the symmetry of the Hamiltonian. For systems without time-reversal 
invariance the relevant ensemble is the Gaussian Unitary Ensemble (GUE), for 
systems invariant with respect to time-reversal the ensemble is the Gaussian 
Orthogonal Ensemble (GOE) and for systems with time-reversal invariance 
but with half-integer spin energy levels have to be described according to the 
Gaussian Symplectic Ensemble (GSE) of random matrices. 

For these classical ensembles all correlation functions which determines 
statistical properties of eigenvalues £"„ can be written explicitly (see e.g. I4(>l . 

The simplest of them is the one-point correlation function or the mean 
level density, d{E), which is the probability density of finding a level in the 
interval {E, E + dE). When d{E) is known one can construct a new sequence 
of levels, e„, called unfolded spectrum as follows 



This artificially constructed sequence has automatically unit local mean den- 
sity which signifies that the mean level density (provided it is a smooth func- 
tion of E) plays a minor role in describing statistical properties of a spectrum 
at small intervals. 

The two-point correlation function, R2{e), is the probability density of 
finding two levels separated by a distance in the interval (e, e + de). The 
characteristic properties of the above ensembles is the phenomenon of level 
repulsion which manifest itself in the vanishing of the two-point correlation 
function at small values of argument 



where the parameter /3 = 1, 2, and 4 for, respectively, GOE, GUE, and GSE. 
This behaviour is in contrast with the case of the Poisson statistics of inde- 
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pendent random variables where 

i?2(e)^ J(i?)^0. 

For later use we present the explicit form of the two-point correlation function 
for GUE with mean density d 

i?2(e) = d' + dS{e) + i?2(e) + 4°"' (e) (40) 
where the smooth part of the connected two-point correlation function is given 

by 

i?2(6) = (41) 

and its oscillating part is 

. (42) 

The term d6{e) in H4UI) corresponds to taking into account two identical levels 
and it is universal for all systems without spectral degeneracy. It is a matter 
of convenience to include it to i?2(e) or not. When one adopts the definition 
(|45() the appearance of such terms is inevitable. 

Another useful quantity is the two-point correlation form factor defined as 
the Fourier transform of the two-point correlation function (unfolded to the 
unit density) 

/C30 
i?2(a;)e2'^**^da; . (43) 
-oo 

For convenience one introduces a factor 27r in the definition of time. 

In Fig. the two-point correlation form factors for usual random matrix 
ensembles are presented. Their explicit formulas can be found in ^Hl- For 
these classical ensembles small-t behaviour of the form factors is 

Kit) (44) 

with the same /3 as above. 

The nearest- neighbor distribution, p{s), is defined as the probability den- 
sity of finding two levels separated by distance s but, contrary to the two-point 
correlation function, no levels inside this interval are allowed. For classical en- 
sembles the nearest-neighbor distributions can be expressed through solutions 
of certain integral equations and numerically they are close to the Wigner sur- 
mise (see e.g. [T^ ^ 

p{s) = as'^e-^^' 

where /3 is the same as above and constants a and b are determined from 
normalization conditions 
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Fig. 8. Two point correlation form factor of classical random matrix ensembles. 

1.5 , ' , ' ^ ' 




''^ • ' • T ^= 1 

12 3 

s 

Fig. 9. Nearest-neighbor distribution for the standard random matrix ensembles. 
Dotted line - the Poisson prediction 



/ p(s)ds — \ sp(s)ds — 1 . 
Jo "'0 

These functions are presented at Fig. [Altogether with the Poisson prediction 
for this quantity p(s) = e~^. 

Though random matrix ensembles were first introduced to describe spec- 
tral statistics of heavy nuclei later it was understood that the same conjec- 
tures can be appUed also for simple dynamical systems and to-day standard 
accepted conjectures are the following 
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• The energy levels of classically integrable systems on the scale of the mean 
level density behave as independent random variables and their distribu- 
tion is close to the Poisson distribution |10) . 

• The energy levels of classically chaotic systems are not independent but 
on the scale of the mean level density they are distributed as eigenvalues 
of random matrix ensembles depending only on symmetry properties of 
the system considered \IJ}. 

- For systems without time-reversal invariance the distribution of en- 
ergy levels should be close to the distribution of the Gaussian Unitary 
Ensemble (GUE) characterized by quadratic level repulsion. 

- For systems with time-reversal invariance the corresponding distribu- 
tion should be close to that of the Gaussian Orthogonal Ensemble 
(GOE) with linear level repulsion. 

- For systems with time-reversal invariance but with half-integer spin 
energy levels should be described according to the Gaussian Symplectic 
Ensemble (GSE) of random matrices with quartic level repulsion. 

These conjectures are well confirmed by numerical calculations. 

The purpose of this Chapter is to investigate methods which permit to ob- 
tain spectral statistics analytically. For a large part of the Section we follow 
|^5|. In Sect. ^ a formal expression is obtained which relates correlation func- 
tions with products of trace formulas. In Sect. II. li the simplest approximation 
to compute such products is discussed. It is called the diagonal approxima- 
tion and it consists of taking into account only terms with exactly the same 
actions. Unfortunately, for chaotic systems this approximation can be used, 
strictly speaking, only for very small time estimated in Sect. 11.21 To under- 
stand the behaviour of the correlation functions for longer time more com- 
plicated methods of calculation of non-diagonal terms have to be developed. 
In Sect. 13 this goal is achieved for the Riemann zeta function. To obtain the 
information about correlations of prime pairs we use the Hardy-Littlewood 
conjecture which is reviewed in Sect. 12.11 The explicit form of the two-point 
correlation function for the Riemann zeros is obtained in Sec. 12.21 In Sect. 13 it 
is demonstrated that the obtained expression very well agrees with numerical 
calculations of spectral statistics for Riemann zeros. 

1 Correlation Functions 

Formally n-point correlation functions of energy levels are defined as the prob- 
ability density of having n energy levels at given positions. Because the density 
of states, d{E), is the probability density of finding one level at point E, cor- 
relation functions are connected to the density of states as follows 

i?n(ei, £2, ...,£„) = {d{E + ei)d{E + £3) ■ ■ • d{E + e„)) • (45) 
The brackets (. . .) denote a smoothing over an appropriate energy window 
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ifiE)) = J f{E')a{E - E')dE' 



with a certain function ij(E). Such smoothing means that one considers eigen- 
values of quantum dynamical systems at different intervals of energy as form- 
ing a statistical ensemble. 

The function (t(E) is assumed to fulfill the normalization condition 



J a{E)dE = 1 



and to be centered around zero with a width AE obeying inequalities 

AEq <t:AE^ AEci < E . 



(46) 



Here AEq has to be of the order of the mean level spacing, AEg « 1/d, and 
AEci denotes the energy scale at which classical dynamics changes noticeably. 
A schematic picture of a{E) is represented at Fig. IIUI 




Fig. 10. Schematic form of smoothing function. 



The trace formula for the density of states of chaotic systems was discussed 
in Chap.^E^nd it has the form 



d{E) = d{E) + J2 \,ne' 



nSp{E)/h 



where the summation is performed over all primitive periodic orbits and its 
repetitions, and 



27r?i|det(Mp"- l)|i/2 



-7rin/ip/2 



(47) 



Substituting this expression in the formula for the two-point correlation func- 
tion one gets 



50 Eugene Bogomolny 
i?2(ei,e2) + 

+ X! ^pi,ni^P2,n2 ( '^'^P Ti^iSp, (E + ei)- naS'pa (E + 62)) ) + c.c. 



Pi ,ni 



and the terms with the sum of actions are assumed to be washed out by the 
smoothing procedure. 

Expanding the actions and taking into account that dS{E) / dE = T{E) 
where T{E) is the classical period of motion one finds 

4^)(ei,e2) = J2 Ap,,n,A;^,„/exp^{n^SpAE)-n2Sp,{E)) 
i 

n 



X exp -(niTpj(£;)ei - n2Tp^{E)e2) + c.c. . 



Here ii2'^^(ei, £2) is the connected part of the two-point correlation function 

i?2(ei,e2) =d' + 4'^(ei,e2). 

The most difficult part is the computation of the mean value of terms with 
the difference of actions 

exp ^(niS'p^(£;) - n2Sp^{E)) 



1.1 Diagonal Approximation 

Berry JI] proposed to estimate such sums in an approximation (called the 
diagonal approximation) by taking into account only terms with exactly the 
same actions having in mind that terms with different values of actions will 
be small after the smoothing. 

Let g be the mean multiplicity of periodic orbit actions. Then the con- 
nected part of the two-point correlation function in the diagonal approxima- 
tion is 

Rf''a\e) ^9 I^P-|'e™^-(^)'/'' + c.c. . (48) 

p. n> 1 

Here e = ei — £2 and the sum is taken over all primitive periodic orbits. 
From (|48|l it follows that the two-point correlation form factor 

/+00 
i?2(e)e2^"Me. 
-00 

in the diagonal approximation equals the following sum over classical periodic 
orbits 

K('^-«)(t) = 27:gJ2\Ap,r.\'S (27:1-^^^^') + c.c. . (49) 
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According to the Hannay-Ozorio de Almeida sum rule sums over periodic 
orbits of a chaotic systems can be calculated by using the local density of 
periodic orbits related with the monodromy matrix, Mp, as follows 

dPp = ^|det(A/p-l)| . 

Using l|T7|) one gets 



gt 



where g is the mean multiplicity of periodic orbits (i.e. the mean proportion 
of periodic orbits with exactly the same action) . For generic systems without 
time-reversal invariance there is no reasons for equality of actions for different 
periodic orbits and 5 = 1 but for systems with time-reversal invariance each 
orbit can be traversed in two directions therefore in general for such systems 
g = 2. Comparing these expressions one concludes that the diagonal approx- 
imation reproduces the correct small-t behavior of form-factors of classical 
ensembles (cf. (I44f) '). 

Unfortunately, K'^'^'^a^t) grows with increasing of t but the exact form- 
factor for systems without spectral degeneracy should tends to d for large t. 
This is a consequence of the following arguments. According to H45|l 

\m,n I 

= S{E - E„)Sie - Em + En)^ . 

If there is no levels with exactly the same energy the second i5-function in the 
right hand side of this equation tends to (5(e) when e — » and the first one 
gives d. Therefore 

i?2(e) ^ d5(e) , whene^O 
which is equivalent to the following asymptotics of the form factor 

K{t) — > d , when i — > cxd . 

This evident contradiction clearly indicates that the diagonal approximation 
for chaotic systems cannot be correct for all values of t and more complicated 
tools are needed to obtain the full form factor. 



1.2 Criterion of Applicability of Diagonal Approximation 



One can give a (pessimistic) estimate till what time the diagonal approxima- 
tion can be valid by the following method. The main ingredient of the diagonal 



52 Eugene Bogomolny 



approximation is the assumption that after smoothing all off-diagonal terms 
give negligible contribution. This condition is almost the same as the condi- 
tion of the absence of quantum interference. But it is known that the quantum 
interference is not important for times smaller than the Ehrenfest time which 
is of the order of 

iE« — ln(l/n), 

where Aq is a (classical) constant of the order of the Lyapunov exponent 
defined in such a way that the mean splitting of two nearby trajectories at 
time t grows as exp(Aoi). For billiards {ka)~^, where a is of the order of 
system size, plays the role of h and Aq = kX where k is the momentum and A 
determines the deviation of two trajectories with length L = kt. The constant 
A which we also called the Lyapunov exponent is independent on k for billiards 
and ^ 

~ In(fco) . 

In the scmiclassical limit fc ^ oo the Ehrenfest time and, consequently, the 
time during which one can use the diagonal approximation tends to zero as 
Infc/fc. 

More careful argumentation can be done as follows. The off-diagonal terms 
can be neglected if 



< 1 . 



But this quantity is small provided the difference of periods of two orbits AT = 
Tpi — Tp.^ times the energy window AE used in the definition of smoothing 
procedure is large 

^(Tpi-TpJZi£»l. (50) 

For billiards Tp = Lp/k and this condition means that one has to consider all 
periodic orbits such that their difference of lengths is 

hk 

Lpi -Lp,:^ — . 

But the number of periodic orbits with the length L for chaotic systems grows 
exponentially 

pAL 

N{Lp <L) = — 

where A is a constant of the order of the Lyapunov exponent. Therefore in the 
interval L,L + SI there is e^^5l/L orbits and the mean difference of lengths 
between orbits with the lengths less than L is of the order of 



AL = Lexp(-AL) 
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To fulfilled the above condition one has to restrict the maximum length of 
periodic orbits, L„i, by 

L„iexp(-Ai„) « . 
In the limit of large with logarithmic accuracy this relation gives 

which corresponds to the maximal time till the diagonal approximation can 
be applied 

™ fc Xk Xk 

As AE <$:E ^ fc2, t,n < tE- 

Another important time scale for bounded quantum systems is called the 
Heisenberg time, tn- It is the time during which one can see the discreteness 
of the spectrum 

tH = 2TTd . 

As for billiards d is a constant 



For the Riemann zeta function the situation is better because (i) in this case 
'momentum' plays the role of 'energy' (the 'action' Ehip is linear in E and 
not quadratic as for dynamical systems) and (ii) the density of states for the 
Riemann zeta function is (ln(i?/27r))/(27r). 
The analog of (|50|l in this case is 

(Inpi - \np2)AE > 1 . 

It means that to apply the diagonal approximation prime numbers have to be 
such that the difference between any two of them obeys 

^-Pae » 1. 

p 

The difference between primes near p is of the order of In p. Hence from the 
above inequalities it follows that diagonal approximation can be used till time 
tm = Inpm where Pm is such that 

Pm ~ AE ■ 

Or with logarithmic precision p,„ < AE. As AE < E (see ()46|l'). pm ~ E and 
the maximum time 

tm ^InE = 2Txd{E) 
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i.e. the diagonal approximation for the Riemann zeta function is valid till the 
Heisenberg time which agrees with the Montgomery theorem |45j . 

This type of estimates clearly indicates that the diagonal approximation 
for chaotic dynamical systems can not, strictly speaking, be used to obtain 
an information about the form-factor for large value of t. Only the short-time 
behaviour of correlation functions can be calculated by this method. (Notice 
that for GUE systems the diagonal approximation gives the expected answer 
till the Heisenberg time but it just signifies that one has to find special reasons 
why all other terms cancel.) 

2 Beyond the Diagonal Approximation 

The simplest and the most natural way of semi-classical computation of the 
two-point correlation functions is to find a method of calculating off-diagonal 
terms. We shall discuss here this type of computation on the example of 
the Riemann zeta function where much more information then for dynamical 
systems is available (for the latter see |22 and '25'). 

The trace formula for the Riemann zeta function may be rewritten in the 
form 



The connected two-point correlation function of the Riemann zeros, i?2 = 
R2 — (P, is 



The diagonal approximation corresponds to taking into account terms with 
ni = 712 




where 






This expression may be transform as follows (cf. |2]) 



,(diag) 
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where 

A{e) = |C(l + ie)P<p(*°''ne) , 

and function <p(d^ag) 

is given by a convergent sum over prime numbers 
.Z>(''»^)(e)^exp('2 5:f;i^cos(mdnp)] . 

\ p m=l / 

In the hmit e 0, C(l + ie) ^ (ie)"^ and ^(*°9)(e) const. Therefore in 
this hmit 

which agrees with the smooth part of the GUE result (|4HI . 
The off-diagonal contribution takes the form 

lE ln(ni/n2)+i(€i In ni— €2 In 712) 

711^712 

The term cxp(i£' ln(rii/rt2)) oscillates quickly if rii is not close to n2- Denoting 

and expanding all smooth functions on r one gets 

^ ^ V + ( e'-/"+-"" ) + c.c. 



47^2 n \ / 



where e — ti — t2- 

The main problem is clearly seen here. The function 

F{n,r) = A{n)A{n + r) 

is quite a wild function as it is nonzero only when both n and n + r are powers 
of prime numbers. As we have assumed that r <C n, the dominant contribution 
to the two-point correlation function will come from the mean value of this 
function over all n, i.e. one has to substitute into r'"2^^^ (e) instead of F(n, r) 
its mean value 



1 ^ 

i{r) = \m\ — A{n)A{n -f r) 



2.1 The Hardy-Littlev^rood Conjecture 

Fortunately the explicit expression for this function comes from the famous 
Hardy-Littlewood conjecture. There are two different methods which permit 
to 'find' this conjecture. We start with the original Hardy-Littlewood deriva- 
tion [55] . 
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First, let us remind two known facts. The number of prime numbers less 
that a given number N{p < x) is asymptotically (see e.g. ^Sj) 

N{p<x) = -^ . 

inx 

Conveniently it can also be expressed in the following form 

N 



lim 

n—1 

The number of prime number Nq^rip < x) in arithmetic progression of the 
form mq + r with (r, q) = 1 and r < g is given by the following asymptotic 
formula (see e.g. PI]) 

X 

NqAP <X)^ —r— 

ip{q) mx 

where i^{n) is the Euler function which counts integers less than n and co- 
prime with n 

(p{n) =nT[ ( 1 - - 

pin- 
As above, this relation can be rewritten in the equivalent form 



N^oo N ^ (p[q) 



m — l 

In the Hardy-Littlewood method one introduces the function 



f{x) = Y,A{n)x- 



n=l 



which converges for all complex x such that < 1. 

In the circle method of Hardy and Littlewood one considers the be- 
haviour of this function close to the unit circle when the phase of x is near a 
rational number 2'iTp/q with co-prime integers p and q. One gets 



\np / q-\-'ui5 



with M, 5 ^ 0. 

In the exponent there is a quickly changing function 2T:np/q. It is quite 
natural to consider n from the arithmetic progression 

n = mq -\- r 

with fixed q and r < g. In this case 
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m. r 

Substituting instead of A{mq + r) its mean value H52|l one gets 



1 f°° 
Vil) , Jo 



In the last step we use that fact that 

where ii{q) is the Mobius function defined through the factorization of q on 
prime factors 

r 1 if 9 = 1 

= < (-1)'' if <? = Pi...Pfe 
[ if g is divisible on 

The final expression means that function f{x) has a pole singularity at the 
unit circle at every rational point. 

The knowledge of f{x) permits formally to compute the mean value of the 
product of two yl-functions. 

Let 

J^{R) = _ / /(i?e''^)/(i?e-''^)e-''''^dv3 = R^'Y + r)A{m)R^"' . 

As the function /(x) has a pole singularity at the unit circle at every rational 
point one can try to approximate this integral by the sum over singularities 



JO 



^r(e-") = ^ / /(i?e"^)/(i?e-'^)e-"-'^d^ 



}_ f j^g-«+i27rp/g+i(5^ j^g-ti-27rip/g-ii5-j^-ir(27rp/ij+ii5) . 

(P, = 1 



27r ^ ^ V ip(q) J .1 u'^ + S'^ 

J_ y^ g27ripr/g f ^ 
(p, q) = l ^^^^^ 



Therefore 



2 



(P, ■/) = ! 
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from which it foUows that 

N 



1 ^ 

hm — A{n)A{n + r) = a{r) 



where 

"M=f:(^)' E (53) 

9=1 ^^'■^''^ (p, g) = l 

Using properties of such singular series one can prove |35| that for even r 
a(r) = and for odd r it can be represented as the foUowing product over 
prime numbers 

"^=^^11^ (54) 

p\r 

where the product is taken over all prime divisors of r bigger than 2 and C2 
is the so-called twin prime constant 

= 211(1- 7^^) « 1-32032.... (55) 

Instead of demonstration the formal equivalence of (|53|l and (|54|l we present 
another heuristic 'derivation' based on the probabilistic interpretation of 
prime numbers which gives directly H54() and (|55|l . 
The argumentation consists on the following steps. 

• Probability that a given number is divisible by a prime p is 

lim — [number of integers divisible hy p < N] — — . 

iV-»oo N p 

In general to find such probabilities it is necessary to consider only the 
residues modulo p and find how many of them obey the requirement. 

• Probability that a given number is not divisible by a prime p is 

1 

1 - - . 

P 

• Probability that a number is not divisible by primes pi,P2, ■ ■ ■ Pk is 

k 



n(i-^) 



(56) 



The above formula is correct for any finite collection of primes but for com- 
putations with infinite number of primes it may be wrong. 
For example, when used naively it gives that 
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• probability that a number a; is a prime is 

lid--)- 

p<y/x 

This prime number 'theorem' is false because from it it follows that the number 
of primes less than x is 



n{x) ^x n (1 ~ -) 7^2e-T 



X — ^oo X 

p' In a; 



which differs from the true prime number theorem by a factor 2e^'^ w 1.123 
where 7 is the Euler constant. The origin of this discrepancy is related with the 
approximation frequently used above: [x/p] = xjp where \x\ is the integer part 
of X. Instead of (|56() one should have Hpll ^ [^^/p]/^)- For a finite number of 
primes and x — > 00 it tends to (j56|l . But when the number of primes considered 
increases with x errors are accumulated giving a constant factor. 

Nevertheless one could try to use probabilistic arguments by forming ar- 
tificially convergent quantities. One has 

• Probability that x and x -V r are primes is 

lim — [ number of integers x < N such that x and x + r are primes 1 . 

N^oo N f I 

Let consider a prime p. Two cases are possible. Either p \r or p/r. In the first 
case the probability that both number x and x + r are not divisible by p is 
the same as the probability that only number x is not divisible by p which is 



p\r 

When pj[r one has to remove two numbers from the set of residues as a; = 
0, 1, ..,p—l (modp) and x+r = 0, 1, ..,p—l (modp). Therefore the probability 
that both numbers x and x + r are not divisible by a prime p is 

Finally 

• Probability that both x and x + r are primes = Y[p\r Hp/r • 

To find a convergent expression we divide both sides by the probability that 
numbers x and x + r are independently prime numbers computed also in the 
probabilistic approximation. The latter quantity is 

[ Probability that x is prime and a; < ] Y\ ■ 
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Therefore 

[ Probabihty that both x and x + r are primes with x,x + r < N ] 
[ Probabihty that x is prime Y 

p\r ph P ^ p>2 ^ ^ p\r 

As the denominator in the above expression is 1/ In^ iV it follows that the 
probability that both a; and x + r are primes with x < N, and a; + r < is 
asymptotically 

a{r) 

with the same function a(r) as in H54() . 

We stress that the Hardy-Littlewood conjecture is still not proved. Even 
the existence of infinite number of twin primes (primes separated by 2) is not 
yet proved while the Hardy-Littlewood conjecture states that their density is 
Ca/ln^iV. 



2.2 Two-Point Correlation Function of Riemann Zeros 

Taking the above expression of the Hardy-Littlewood conjecture as granted 
we get 

Ri°"\e) = ^ E ^e"'"" E «(r)e'^'-/" + c.c. . 

n>l r 

After substitution the formula for a{r) and performing the sum over all r one 
obtains 

where the summation is taken over all pairs of mutually co-prime positive 
integers p and q (without the restriction p < q). 

Changing the summation over n to the integration permits to transform 
this expression to contributions of values of n where 

^-^ = 0. 
q zirn 

In this approximation 




+ c.c. . 



Using the formula (which is a mathematical expression of the inclusion- 
exclusion principle) 
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oo 

(p,q) = l k=l S\q 

and taking into account that 27ro? — ln(i?/27r) one obtains 

R^°ff^ (e) = ^ 1^(1 + i£) |2e2-irf^<^(o//) (e) + c.c. (57) 
where function ^(°-^-^) (e) is given by a convergent product over primes 

and <p(°//)(0) = 1. 

In the Umit of small e 

which exactly corresponds to the GUE results for the oscillating part of the 
two-point correlation function H42|l . 

The above calculations demonstrate how one can compute the two-point 
correlation function through the knowledge of correlation function of periodic 
orbit pairs. For the Riemann case one can prove under the same conjectures 
that all n-point correlation functions of Riemann zeros tend to corresponding 
GUE results 

3 Summary 

Trace formulas can formally be used to calculate spectral correlation functions 
for dynamical systems. In particular, the two-point correlation function is the 
product of two densities of states 

i?2(e) = {d{E + e)d{E)) . 

The diagonal approximation consists of taking into account in such products 
only terms with exactly the same action. For chaotic systems this approxima- 
tion is valid only for very small time. In particular, it permits to obtain the 
short-time behaviour of correlation form factors which agrees with predictions 
of standard random matrix ensembles. 

The main difficulty in such approach to spectral statistics is the necessity to 
compute contributions from non-diagonal terms which requires the knowledge 
of correlation functions of periodic orbits with nearby actions. 

For the Riemann zeta function zeros it can be done using the Hardy- 
Littlewood conjecture which claims that the number of prime pairs p and 
p -\- r such that p < N for large N is asymptotically 
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air) — 5 — 

where a(r) (with even r) is given by the product over aU odd prime divisors 
of r 

and 

p>2 ^ ' 

Using this formula one gets that the two-point correlation function of Riemann 
zeros is 



where the diagonal part 



1 d 



2 



and non-diagonal part 

R[°^f\e) = ^IC(1 + ie)|'e2"''^"'<?(°^/)(e) + c.c. . 

The functions ^>(''^"ff)(e) and ^"ff ^e) are given by convergent products over 
all primes 

^(*-)(e) =exp £ i^cos(mdnp)) 

\ p m=l ^ / 

and 

In a few other methods were developed to 'obtain' the two-point correla- 
tion function for Riemann zeros. These methods were based on different ideas 
and certain of them can be generalized for dynamical systems. Though neither 
of the methods can be considered as a strict mathematical proof, all lead to 
the same expression (I57II . 

It is also of interest to check numerically the above formulas. When numer- 
ical calculations are performed one considers usually correlation functions for 
the unfolded spectrum. For the two-point correlation function this procedure 
corresponds to the following transformation 

^(unfolded) = 3^ ^2 



d^E) - \d{E) 
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At Fig II II we present the two-point correlation function for 2 • 10* zeros near 
the 10^^-th zero computed numerically by Odlyzko 0H| together with the GUE 
prediction for this quantity 

At Figs. I12I15I we present the difference between the two-point correlation 




12 3 



Fig. 11. Two point correlation function of the Riemann zeros near the 10^3-th zero 
(dots) and the GUE prediction (solid line). 

function computed numerically and the GUE prediction. At Fig. E| we 
present the difference between numerically computed two-point correlation 
function and the 'exact' function and at Fig. El the histogram of differences 
is given. Notice that these differences are structure less and the histogram 
corresponds practically exactly to statistical errors inherent in the calculation 
of the two-point correlation functions which signifies that the obtained formula 
agrees very well with the numerics. 
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0.004 I ^ 




Fig. 12. The difference between the two point correlation function of the Riemann 
zeros and the GUE prediction in the interval < e < 5. The solid line is the 
difference between the 'exact' correlation function and the GUE prediction in this 
interval. 
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Fig. 14. The same as at Fig. 1121 but in tlie interval 10 < e < 

0.004 I ^ 1 ^ 1 ^ 1 ^ 1 ^ 1 
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Fig. 15. The same as at Fig. 1121 but in the interval 15 < e < 
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e 

Fig. 16. The difference between numerically computed two-point correlation func- 
tion and the 'exact' function 
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Fig. 17. The histogram of the deviations of the numerically computed two-point 
correlation function of Ricmann zeros and the 'exact' formula. Solid line is the 
Gaussian fit to the histogram. 
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Arithmetic Systems 



As was discussed above it is well accepted that spectral statistics of classically 
chaotic systems in the universal limit coincides with spectral statistics of the 
usual random matrix ensembles. But it is also known (see e.g. |7], |29| and 
references therein) that the motion on constant negative curvature surfaces 
generated by discrete groups (considered in Chap. ^ is the best example of 
classical chaos. Consequently, models on constant negative curvature seem to 
be ideal tools to check the conjecture on spectral fluctuations of classically 
chaotic systems. Their classical motion is extremely chaotic and time-reversal 
invariant and a priori assumption was that all of them should have energy 
levels distributions close to predictions of the Gaussian orthogonal ensemble 
(GOE) of random matrices. 

Nevertheless when the first large scale numerical calculations were per- 
formed [3], [H^ they clearly indicated that for certain hyperbolic models the 
spectral statistics were quite close to Poisson statistics typical for integrable 
systems. 

As an example we present in Fig. llSl the nearest-neighbor distribution for 
the hyperbolic triangle with angles (7r/2, 7r/3, 0) corresponding to the well- 
known modular triangle with Dirichlet boundary conditions. The agreement 
with Poisson prediction is striking though classical motion for this system is 
perfectly chaotic. 

The purpose of this Chapter is to show that this strange behaviour is the 
consequence of exponentially large exact degeneracy of periodic orbit lengths 
in systems considered 18^. In all hyperbolic surfaces where the Poisson-like 
statistics was observed there is on average exp(Z/2) classical periodic orbits 
with exactly the same length I. It will be demonstrated that this is the char- 
acteristic property of models generated by the so-called arithmetic groups. 
As classical mechanics is not sensitive to lengths of periodic orbits all these 
models remain completely chaotic. But the cumulative effect of interference 
of degenerated periodic orbits changes drastically the quantum mechanical 
properties. 
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Fig. 18. The nearest neighbor distribution for 10000 first levels of the triangle 
(7r/2, 7r/3, 0) (the modular triangle). Solid line - the Poisson distribution. Dotted 
line - the GOE distribution. 

This Chapter is based mostly on 24 . In Sect. simple calculations prove 
exponential degeneracy of periodic orbit lengths for the modular group. The 
main peculiarity of the modular group matrices is that their traces are inte- 
gers. Therefore if one considers all matrices with |Tr M\ < X the number of 
different traces increases at most linearly with X. In Sect.|21it is shown that 
this property is typical for all arithmetic groups. An informal mini-review of 
such groups is given in this Section and it is demonstrated that for all these 
groups exponentially many periodic orbits have exactly the same length. From 
the results of it follows that there is exactly 85 triangles generated by dis- 
crete arithmetic groups. All triangular models were the Poisson- like spectral 
statistics was numerically observed are in this list. In Sect.|31it is shown that 
in the diagonal approximation the two-point correlation form factor of arith- 
metic systems jumps very quickly to the Poisson value thus confirming unusual 
nature of arithmetic systems. In Sect.01the exact two-point correlation func- 
tion for the modular domain is calculated. The correlations of multiplicities 
are obtained by a generalization of the Hardy-Littlewood method discussed 
in Sect. 12. ll The resulting formula proves that in the universal limit the two- 
point correlation function of eigenvalues of the Laplace-Beltrami operator au- 
tomorphic with respect to the modular group tends to the Poisson prediction. 
Arithmetic groups have many other interesting properties. In particular, for 
all arithmetic groups it is possible to construct an infinite number of mutually 
commuting operators which commute also with the Laplace-Beltrami opera- 
tor. Properties of these operators called the Hecke operators are discussed in 
Sect.|5| The Jacquet-Langlands correspondence between different arithmetic 
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groups is mentioned in Sect. In Sect. [7| non- arithmetic models are briefly 
discussed. 



1 Modular group 

The modular group is the group of all 2 x 2 matrices 

M ■ 



m n 
k I 

with integer m, n, k, I and the unit determinant ml — nk = 1. 

The periodic orbits correspond in a unique way to the conjugacy classes of 
hyperbolic elements of the group (see Sect. I2.8|l . The length of periodic orbit 
Ip is related with the trace of a representative matrix of the conjugacy class 
M as follows 

2 cosh ^ = iTr Ml . 
2 ' ' 

As all elements of modular group matrices are integers, the trace is also an 
integer 

|TrM| = n . (58) 

Here the arithmetical nature of the group clearly comes into the play. This 
simple property is very important. It signifies that for the modular group there 
is just a quite restrictive set of all possible traces and, consequently, of periodic 
orbit lengths. For modular group the number of possible different lengths is 
the number of different integers less than 2coshL/2 (see (jSHl)): hence 

Ndif. lengths = 2 cosh — ye'. 

On the other hand, for any discrete group the number of periodic orbits of 
length less than L grows asymptotically as 

Let g{l) be the multiplicity of periodic orbits with length I. One has obvious 
relations valid for large L 

KL KL 

where the summation extends over different lengths of periodic orbits counted 
without taking multiplicity into account. 

Let us define the mean multiplicity {g{l)) as the following ratio 

Number of periodic orbits with I < Ip < I + Al 
[9( )) - jv^y^i^gj. different lengths with I <lp <l + Al' ^ ' 
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Asymptotically for large L the previous formulas gives 

(.9(0) = 2— 

which demonstrates that periodic orbit lengths for the modular group are 
exponentially degenerated. 

2 Arithmetic Groups 

The crucial feature which led to the exponential degeneracy of periodic orbit 
lengths for the modular group was the fact that traces of modular group 
matrices were integers which was a direct consequence of the arithmetic nature 
of modular group. But 2x2 matrix groups with integer elements are exhausted 
by the modular group and its subgroups. 

Nevertheless, one can construct a quite large class of discrete groups with 
strong arithmetic properties by considering groups which are not equal to 2 x 2 
integer matrices but which permit a representation by n x n integer matrices 
(n > 2). 

The existence of such representation means that for each 2x2 group 
matrix, g, one can associate a n x n matrix with integer entries, M{g), in 
such a way that the matrix associated to the product of two groups matrices 
equals the product of two matrices associated to the corresponding factors 

M{ab) = M{a) x M{b) 

for all a and b from the group considered and M{1) = 1. 

To define general arithmetic groups we need a few definitions. 

• A subset of a group F is called a subgroup if it forms itself a group. 

• A subgroup g of a group F is called a subgroup of finite index if F can be 
represented as a finite union 

r ^ g + g^i + ■ ■ ■ + gjk 

with 7fc e P. 

• Two groups are called commensurable if they have a common subgroup 
which is of finite index in both of them. 

Groups which have a representation by integer matrices and all groups 
commensurable with them are called arithmetic groups. This Section is de- 
voted to the investigation of their properties. 

In Sect. l2.l1 a non- formal review of algebraic fields is given and in Sect.O 
the construction of quaternion algebras over algebraic fields is shortly dis- 
cussed. It appears that all arithmetic groups can be obtained from quaternion 
algebras with division and in Sect. the necessary and sufficient conditions 
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that a given group will be an arithmetic group is presented. Using these con- 
ditions in Sect. 12.41 it is proved that periodic orbit lengths for all arithmetic 
groups have the same exponential degeneracy (up to a constant factor) as for 
the modular group. 

2.1 Algebraic Fields 

Everybody is familiar with usual rational numbers 

P 

u — — 

q 

with integer p and q. Their important properties are (i) that the sum and 
the product of any two rational numbers also have the same form and (ii) all 
elements except have an inverse (i.e. the division is always possible). From 
mathematical viewpoint rational numbers form a field called Q. 

Algebraic fields of finite degree, IF, are a generalization of this reference 
field obtaining by adding to the set of rational numbers a root a of an irre- 
ducible polynomial 

n 

J2 cfe"'' = (60) 

fe=0 

with integer coefficients c^. This field is denoted IF = Q(a). 
Each element u €;Q(q;) can be represented by the sum 

1=0 

where the bi are usual rationals. The summation and the multiplication of 
these elements are done as with usual numbers except that all powers of a 
larger than n — 1 have to be reduced using the defining equation H6U|) . 

Integers of the field Q(a) are its elements which obey a polynomial equa- 
tion with integer coefficients with an additional condition that the highest 
power coefficients equals one (such polynomials are called monic polynomi- 
als). 

In general, algebraic integers, u, of a field of degree n are freely generated 
by n linearly independent elements of the field Pk with integer coefficients (in 
mathematical language it means that they form a free ffi- module of rank n). 
Explicitly 

n-l 

w = ^ mkPk (61) 

fc=0 

where all nik are usual integers. In simple cases /3k = and 

u) = mkCt^ 

k=0 
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with integer or demi-integer coefficients m^j. Algebraic integers like usual in- 
tegers form a ring (not a field) because the division is not always possible. 

The polynomial equation defining the field H6U|I has n different roots ai,i = 
0, 1, n — 1 with Q!o = a. Any relation between elements of the field remains 
unchanged under the transformations 



where one substitutes in all expressions instead of one root a another root 
Qj. These transformations are called isomorphisms or embeddings of this field 
into(D and they are the only transformations respecting the laws of the field. 

Example 

Add to the field of rational numbers Q one root of the quadratic equation 

=d 

where d is a square- free integer. Elements of this field Q(\/d) can be written 
as 

u = p + qV2 

with p and q rationals. Let 

Lo = a + bV2 

be integers of this field. To find values of a and b one notes that ui obeys the 
quadratic equation 

- 2auj + 0^ - db'^ ^0 . 

To describe integers 2a and a^ — db^ have to be usual integers. Depending on 
d two types of solutions are possible. 

• Ifd=2ord = 3 (mod 4) then a and b have to be integers and 

u! = m + ny/d 

with integers m and n. 

• If d = 1 (mod 4) then both a and b can be demi-integers and 

m n r- 

u) = 1 V d 

2 2 

with integers m = n (mod 2). 

To avoid the last restriction this expression can be rewritten in the form 

l + Vd 

uj = m + n ^ (62) 

with arbitrary integers m and n. 
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As this field is defined by an equation of the second degree it has two isomor- 
phisms 

(j>o ■ p + qVd — > p + qVd , 
(j>i : p + q\/d — > p - qVd . 

Because the product of two algebraic integers is also an algebraic integer from 
H61|l it follows that all algebraic integers permit a representation by matrices 
with integer elements in such a way that the matrix representing a product 
of two integers equals the product of matrices representing each factor. 

For example, for the above considered case of Q(Vrf) with 1 (mod 4) 
one can associate with an integer of this field, oj = m + n\/2, a 2 x 2 matrix 

It is easy to check that this is the true representation of field integers because 
M(wiW2) = M(wi)M(w2) and M(l) = 1. 

When d = 1 (mod 4) the integers have the form (|62() and one can check 
that the matrix representation can be chosen as follows 

2.2 Quaternion Algebras 

Algebras are more general objects then fields. A (vector) algebra of finite 
dimension d is defined as formal sum 

7 = Xil2 + X2h + • • • + XdU ■ 

Here Xj belong to a basis field IF and ij are formal objects (vectors) with a 
prescribed multiplication table 

d 
p=i 

where Cj^. are from the basis field. The sum and the product of any two 
elements of an algebra belong to it. General algebras should be neither com- 
mutative, nor associative. 

An algebra is called a normed algebras if there exists a function, N{'j), 
which associates to any element of the algebra a number from the basis field 
such that the norm of the product equals the product of the norms of both 
factors 

iV(7i72) = N{-fi)N{-f2) . 
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An algebra is called a division algebra if the division is always possible (except 
a zero element). 

Finite dimensional normed division algebras over real numbers are ex- 
hausted by the following three possibilities (the Frobenius theorem j44)'). 

• Commutative and associative division algebras are isomorphed either to 
the usual field of real numbers IR or to the field of complex numbers (D. In 
the latter case the algebra is given by 

7 = .Tl + X2i 

with i^ = —1. The norm in this case is 

N{j) = xl + xl . 

• Non-commutative but associative division algebras are isomorphed to the 
quaternion algebra 

^ ^ xi + X2i + xsi + Xik (64) 

where 

i2=f =k2 = -l, k = ij = -ji. (65) 
The norm of the quaternion algebra is 

iV(7) = xl+ xl + xl + xj . 

• Non-associative normed division algebras are isomorphed to the octonion 
algebra 

8 

fc=l 

with a complicated multiplication table and the norm given by the sum of 
8 squares 



xl 



k ■ 

fe=l 



Similarly, for an algebraic field IF of finite degree there exist quaternion 
normed algebras defined similarly to Eqs. H64(l and 165(1 '56'. These alge- 
bras are labeled by two elements a, 6 e IF and it is a four-dimcnsional non- 
commutative algebra with basis (1, i, j, k) as in ((64(1 with the following multi- 
plication table 

i2=a, f = b, k = ij--ji. (66) 
Such algebra is denoted by {^j^^ and its norm is 

iV(7) = - axl - bxl + abxl . (67) 

The matrix representation of the quaternion algebra 1(66(1 is obtained by the 
isomorphism 
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^ \ . fOl\ , .. / ^^ 

k = ij 



-V^ J ' \bOJ ' \~b,/E 

Explicitly 

b{x3 - Xi^/a) xi - X2Va J 

with X2, 2^3, 2:4 G IF. As it is a representation of the quaternion algebra the 
product of these matrices also has the same form. In this representation the 
norm of the algebra (|67|) equals the determinant of the matrix (|68l) 

N{j) = det 7 . 

From an algebraic field IF of finite degree one can build also another simple 
set of matrices called M(2,IF) given by 2 x 2 matrices with entries from IF 



Xi X2 
X3 Xi 



(69) 



Are the two sets (I68|l and (I69|) different or are they isomorphic? For example, 
if a = and w e F the set is, evidently, within M(2,F). 

Let us show that if ^/a ^ W but there exist certain elements qi, q2, qs, Qi G 
W such that the determinant of the matrix H68(l is zero 

det(7) = qf- qla - h{ql - qla) = (70) 

then matrices H68|) are isomorphic to M(2,IF) |42j . 

Indeed, from the above expression it follows that in this case 6 has the 
form 

b= {ql - qla){q^ - q\aY^ = {ui + U2^/a){ul - U2\/a) 

where 

ui + U2\fa = (qi + q2\fa){(iz + q^-JaY^ . 

As all fractions of field elements belong to the field IF, u\ and U2 are also 
elements of W . Now one can check that 



X\+X2\fa xz-Vx^^fa 
b{x3 - Xiy/a) x\- X2\fa 



where 



Xi + X3U1 + X4U2a X2 ~ X3U2 ~ XaU\ 

a{x2 + X3U2 + XiUi) xi — X3U1 — X4,U2a 



M ^ 

and S' is a fixed (independent of Xi) matrix 
S 



ui+U2\/a 1 



The importance of such representation lies in the fact that the matrix M 
contains only elements of our basis field IF and docs not contain ^/a. Therefore, 
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it belongs to M(2,IF) and the set of matrices 7 (|68|l is the conjugation of 
matrices from M(2,IF) by a fixed matrix S. In other words, when the equation 
det(7) = has a solution 7 S IF expression (|68|l is just a complicated way 
of writing matrices from M(2,IF). These considerations demonstrate that in 
order to construct a group different from M(2,IF) it is necessary to require 
that there exist no elements from the basis field such that the determinant 
()7()|l equals zero. Or, equivalently, any matrix I|fi8(l should have an inverse 
element. In the language of quaternion algebra this property corresponds to 
the division algebra for which any element has an inverse. 

As for real fields this condition is quite restrictive and an explicit answer 
can be obtained only in simple cases. Let us consider for example the case 
when IF is the field of usual rational numbers IF =Q. The following theorem 
1321 gives a series of division algebras overQ. 

Let 6 be a prime number and a be an integer such that the equation 

= a (mod b) 

has no integer solution. Then the pair (a, 6) defines a division algebra overQ 
or equivalently the equation (|7()(l 

^1 ~ ■^2'^ ~ ^(^3 ~ 2:40) — (71) 

has only zero rational solutions. 

To prove the theorem note that due to homogeneity of this equation it is 
sufficient to consider integer solutions xi, X2, x^, X4 without common factors. 
From l(7T|) it follows that 

Xi = X2a (mod b) . 

Consider first the case when b does not divide X2, 6/a;2. As b is assumed to 
be a prime, a;^^ (mod b) exists and {xi/x2)^ = a (mod b) which contradicts 
our assumption. Hence b \ x^ but then b \ x\ and 

x\ = x\a (mod &) . 

The same arguments give that b \ x^ and b \ 2:3 which contradicts the assump- 
tion about the absence of common factors of xi . Therefore there is no rational 
solution of l|71|) and the quaternion algebra defined by a and & is a division 
algebra. 

Quaternion algebras with division are analogs of algebraic fields. How one 
can define integers of a quaternion algebra? 

We have seen above that algebraic integers of a field of degree n form a 
free ffi-module of rank n i.e. they can be represented as a sum of n elements of 
the field with integer coefficients (see (jfilll . Similarly one can define 'integers' 
of a quaternion algebra over such field as a free K-module of rank 4n (which 
generates the whole algebra). For technical reasons they are called not integers 
but 'the order' in the algebra. The word 'integers' in algebras is reserved for 
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elements for which the trace and the determinant of matrix (|68|) are integers 
of the basis field. Different orders exist and the one which is not contained in 
any other order is called the maximal order. 

The simplest case appears when a and b are integers of the basis field IF. 
Then matrices of the form 

xi + X2\/a X3 + XA\fa 
b{x3 — X4^) xi — X2\fa 

where all Xk are integers of IF form an order of the algebra {^^^ (but not 

necessarily the maximal order). 

Matrices of the order in a division quaternion algebra with unit deter- 
minant form a group. Each matrix of this group belongs to the order and, 
therefore, is defined by 4?! integers. The product of two group matrices have 
the same form and corresponds to a certain transformation of integers defin- 
ing both matrices. It means that these groups can be represented by 4n x 4n 
matrices with integer elements. 

All such groups, all their subgroups, and all groups commensurable with 
them are discrete arithmetic groups with finite fundamental domain |42|. 

Example 

As a;^(mod 5) takes only values 0, 1, 4 the equation 

x^ = 3 (mod 5) 

has no integer solution. Hence the pair (3,5) defines a division algebra overQ. 
A simplest order of this algebra has the form 

f m + nVS k + l^/i \ , , 

\h{k - m - nVi ) ^ ' 

with integer m,n,k,l. When one considers these matrices with the unit de- 
terminant 

- Sn^ - 5k^ + IbP = 1 

they form a discrete arithmetic group Fi with a finite fundamental area. 

The order 1)72(1 is not the maximal order. The latter can be chosen e.g. as 
follows 

f ^{m + nVi) Uk + lVS) \ 

with integer to, n, fc, I such that m = k (mod 2) and n = I (mod 2). Matrices 
H73|) with the unit determinant 

- 3n^ - 5fc2 + 151^ = 4 



constitute another discrete arithmetic group /2 whose fundamental domain 
is smaller than for the group (|72|) as Ji is a subgroup of /2. 
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Using the representation Ht)3|l one concludes that to each 2x2 matrix of 
the group Fi one can associate the 4x4 matrix with integer elements 



It is straightforward to check that (i) M(7i72) — M(7i)M(72) for all 
71 , 72 e A, (ii) M(l) = 1, and (iii) det(M) = (det(7))2 = 1 . Together 
these expressions mean that this group is an arithmetic group. 

2.3 Criterion of arithmeticity 

For general fields the situation is more complicated. To explain the general 
criterion of arithmetic groups let us first stress a difference between usual 
integers and algebraic integers. 

The usual integers correspond to a discrete set of points. But for general 
algebraic integers this is not the case. For example, in the field Q(-\/2) integers 
have the form n + m\/2 with integer n and m. But it is evident that one 
can construct sequences of these algebraic integers converging to zero, e.g. 
(\/2 — 1)''' = A/fc — NkV^ when k — > oo. Therefore the set of algebraic 
integers is not discrete as it have finite accumulation points. 

How can one deal with such problem? The main point is that these small 
terms become large under the transformation 



It is clear that now arbitrary small integers are excluded and one gets a discrete 
set of points. 

For more general fields the transformation H74|) is generalized to all non- 
trivial isomorphisms of the field. To remove arbitrary small elements one has 
to require that for all isomorphisms of the field (except the identity), (pi, 
transformed values of integers Ht)l|) are restricted 



M(7) 



^ m n k l\ 
2n m 21 k 
5fc —51 m —n 

\-10k bk-2n nij 



V2 



V2 



(74) 



Let consider in the above example not all algebraic integers n + m\/2 but only 
those which after transformation (|74|l remain bounded 

V2\ < constant . 




(75) 



In order to be sure that all small numbers are removed it is necessary that 
all roots of defining equation H60I) are real. Otherwise, changing a root to its 
complex conjugate may not change modulus of integers. 
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These considerations make reasonable that in order to construct a discrete 
subset of algebraic integers (without finite accumulation points) it is necessary 
that (i) the field be a totally real field (i.e. all roots of defining equation H6(J|I 
are real) and (ii) for all non-trivial isomorphisms of the field transformed 
integers remain bounded as in lf75|l . 

A precise criterion of arithmeticity obtained by Takeuchi |53| is quite sim- 
ilar (see also [H] and [21] for particular examples). 

Takeuchi proved that a group F is an arithmetic group if and only if the 
traces of group matrices have the following properties 

• All Tr(7) are integers of a totally real algebraic field of finite degree. 

• For any non- trivial isomorphism of this field which changes some |Tr(7)| 
for certain 7, the value of the transformed trace satisfies |0(Tr(7))| < 2. 

There are two types of arithmetic groups. Non-compact groups, built from 
SL{2,7L), and compact ones built from quaternion algebra different from 
M(2,Q). 

The above criterion is quite effective, in particular, it permits to find all 
possible arithmetic groups with triangular fundamental domains [53] . There 
are 85 triangular hyperbolic surfaces generated by discrete arithmetic groups. 
All of them are given in Tabled 

2.4 Multiplicities of Periodic Orbits for General Arithmetic 
Groups 

The geometrical length of the periodic orbit, /, is connected with the trace of 
class of conjugate matrices by (I28f) . When / oo 

exp^ = |Tr(7)| . 

Let us prove that for an arithmetic group the number of possible values of 
group matrix traces obeys 

7V(|Tr(7)| <R)^-^ C-R 

with a constant C depending on the group. The traces of matrices of arithmetic 
groups are dispatched as usual integers among real numbers. 

Let r be an arithmetic group. The set of traces {Tr(7), 7 £ F} are integers 
of an algebraic field IF 

1=0 

where rrii are integers and Pi are linearly independent elements of the field. 
Consider the simplest case (3i = a' then 

n-l 

Tr(7) ^to^J2 • 

i=0 



80 Eugene Bogomolny 



Table 1. The list of arithmetic triangles from |53) . {n,m,p) in the first column 
corresponds to the three angles (t^ jn^-K Im,!^ j'p). The second column indicates the 
algebraic field from which is built the corresponding arithmetic group. 



(m,n,p) 


IF 


(2,3,00) 

(3,00,CX3) 


(2,4,cx)) (2,6,00) (2,00,00) 
(4,4,00) (6,6,00) (00,00,00^ 


(3,3,00) 


Q 


(2,4,6) 


(2,6,6) (3,4,4) (3,6,6) 




Q 


(2,3,8) 
(3,8,8) 


(2,4,8) (2,6,8) (2,8,8) 
(4,4,4) (4,6,6) (4,8,8) 


(3,3,4) 


Q(^) 


(2,3,12) 
(6,6,6) 


(2,6,12) (3,3,6) (3,4,12) 


(3,12,12) 


Q(\/3) 


(2,4,12) 


(2,12,12) (4,4,6) (6,12,12) 




Q(V3) 


(2,4,5) 
(5,10,10) 


(2,4,10) (2,5,5) (2,10,10) 


(4,4,5) 


Q(\/5) 


(2,5,6) 


(3,5,5) 




Q(\/5) 


(2,3,10) 


(2,5,10) (3,3,5) (5,5,5) 




(D(a/5) 


(3,4,6) 


(B(\/6) 


(2,3,7) 
(3,3,7) 


(2,3,14) (2,4,7) (2,7,7) 
(7,7,7) 


(2,7,14) 


Q(cos7r/7) 


(2,3,9) 
(9,9,9) 


(2,3,18) (2,9,18) (3,3,9) 


(3,6,18) 


Q(cos7r/9) 


(2,4,18) 


(2,18,18) (4,4,9) (9,18,18) 




Q(cos7r/9) 


(2,3,16) 


(2,8,16) (3,3,8) (4,16,16) 


(8,8,8) 


Q(cos7r/8) 


(2,5,20) 


(5,5,10) 




Q(cos7r/10) 


(2,3,24) 
(12,12,12) 


(2,12,24) (3,3,12) (3,8,24) 


(6,24,24) 


Q(cos7r/12) 


(2,5,30) 


(5,5,15) 




Q(cos7r/15) 


(2,3,30) 


(2,15,30) (3,3,15) (3,10,30) 


(15,15,15) 


Q(cos7r/15) 


(2,5,8) 


(4,5,5) 




Q(V2, ^fl) 


(2,3,11) 


Q(cos7r/ll) 



For a field of degree n there exist n — 1 non-trivial isomorphisms : a ak 
where ak is a root of the defining polynomial different from a. 

Suppose that all such transformations change |Tr(7)|. According to the 
criterion of Takeuchi all transformed traces satisfy 

\tk\ < 2 

where 

n-l 

tk = 0fe(Tr(7)) = ^ TUial . 

i=Q 

Consider these equations as transformations from variables ti to new variables 
rrii 122] ■ The volume elements in these two representations are related as 
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dtodti...dtn-i = I J'|dTOodmi...dm„_i 

where 

J = det 



\ dnik ) 



is called the discriminant of the field and in our case (when jdi — a') 
J = det(a^j^)|fc,j=o,...,n-i = n^"' ~ ■ 

i¥=J 

As nii are integers the volume of the smallest cell is one, and the total number 
of possible integer solutions is asymptotically 

A^(|Tr(7)| <R)^ iV(|to| < i?, < 2) ~ C • i? 

where C = /J . 

For any surface of finite area generated by a discrete group the total num- 
ber of periodic orbits with length less that a given value is asymptotically the 
following 

N,^,{1 < L) ^ . 

The number of periodic orbits with different lengths is the same as the number 
of group matrix traces 

^diff. icngths(/ <L)^C- e^/2 
Let g{l) be the multiplicity of periodic orbits with length I. Then 

^g(0-^and ^ 1 = Ce^^ 

KL KL 

where the summation is done over different lengths. 

Finally the mean multiplicity of arithmetic systems defined as in (|59|l has 
the following asymptotics 

Thus we demonstrate that the arithmetic nature of arithmetic groups leads 
to exponential multiplicities of periodic orbit lengths. 

For generic systems one usually does not expect any degeneracy of pe- 
riodic orbit lengths except the ones which follow from exact symmetries of 
the model. For example, systems with time-reversal invariance, in general, 
should have the mean multiplicity equal to 2, which corresponds to the same 
geometrical periodic orbits spanned in two directions. Therefore, arithmetic 
systems are very exceptional in this respect as they display exponentially large 
multiplicities of periodic orbit lengths. Notice, nevertheless, that according to 
Horowitz- Randol theorem |41| . [Slj this degeneracy is unbounded for any sur- 
face generated by a discrete group. However degeneracies of this theorem are 
much smaller than exponential. 
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3 Diagonal Approximation for Arithmetic Systems 

The large multiplicities of periodic orbit lengths in arithmetical systems seem 
to have no importance in classical mechanics. These systems are as chaotic as 
any other models of free motion on constant negative curvature surfaces with 
finite area. Nevertheless, the quantum spectra of these systems are anomalous: 
is it connected to these degeneracies? In this Section we estimate the quan- 
tum two-point correlation form factor for arithmetic systems in the diagonal 
approximation as was done in Sect. [TTTl for generic chaotic systems. 

Assume that there exist g{l) periodic orbits with the same length I. Exactly 
as it was done in Sect. II. ll one gets the following expression for the diagonal 
approximation of the two-point correlation function 

Rf'''\e) = lAMfgilpy"-^'' + c-c. (77) 

where the summation is done over all periodic orbits. The only difference 
with H48(l is that in Sect. ll.ll it was assumed that g is a constant but here the 
multiplicity g = g{l). 

Define the two-point correlation form factor as the Fourier transform of 
i?2(e) 

+ 00 

ite 



/ + 00 
i?2(e)e^ 
-oo 



This definition differs from the previous one by the absence of the factor 2tt 
in the exponent. For later purposes it is more convenient. 

Equation H77|l leads to the following expression for the two-point correla- 
tion form factor in the diagonal approximation 

^(<i»9)(i) = 2^^|Ap^„(giVg'5('t~^') . (78) 

From H76|l it follows that the mean multiplicity of periodic orbit lengths for 
arithmetic systems is asymptotically 

2e'p/2 

^^^'^^^ = ^7 

with a model dependent constant C (for the modular group C — 1). 

For any models generated by discrete groups the summation over all peri- 
odic orbits is asymptotically equals the integration with the following measure 



E 



dl I 

r 



Taking into account that when / — > oo the term with n = 1 dominates and 
(see (ITl ) 
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Ank 

one obtains that in the diagonal approximation 

kt 

It means that the correlation form factor K{t) for arithmetic systems grows 
much faster than was usually assumed and that for time of order of the Ehren- 
fest time it becomes of the order of 1. 

The simplest approximation to the full form factor is the following 



K{t) 



Kdiag^^^ for t < t* 

d for t > t* 



where t* is defined by the requirement that K'^'^°-3{t*) — d 

t* ^- ln(27rfcCJ) . 
k 

For the true Poisson statistics K{t) always equals d. For usual integrable 
systems K{t) increases to this value during the time of the order of shortest 
periodic orbit periods, t* ^ 1/k. For arithmetic systems K{t) jumps to the 
universal saturation value in a time of order of the Ehrenfest time which has 
an additional logarithm of the momentum. 

Therefore, spectral statistics of arithmetic systems is much closer to the 
Poisson prediction typical for integrable systems than to any of standard ran- 
dom matrix ensembles conjectured for generic ergodic systems. 



4 Exact Two-Point Correlation Function for the 
Modular Group 

The diagonal approximation gives quite crude estimate of the form factor. 
For the modular group it is possible to compute explicitly the two-point cor- 
relation function The calculations are based on a generalization of the 
Hardy-Littlewood method and depend strongly on the number-theoretical 
properties of the multiplicities of the periodic orbits of the modular group. 
In Sect. 14. l1 using the Selberg trace formula the two-point correlation form 
factor is expressed through the two-point correlation function of multiplicities 
of periodic orbit lengths for the modular group. In Sect. 14.21 the latter is cal- 
culated by a certain generalization of the Hardy-Littlewood method. Quite 
tedious explicit formulas are given in Sect. 14.31 and the final expression for the 
two-point correlation form factor is presented in Sect. ^31 
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4.1 Basic Identities 

The modular group has been considered in Sect-Q] It is the group of all 2 x 2 
matrices with integer elements and unit determinant. The periodic orbits of 
the modular group correspond in a unique way to the conjugacy classes of 
hyperbolic elements of the modular group. The length of periodic orbit Ip is 
related with the trace of a representative matrix of the conjugacy class as 
follows 

|TrAf| = 2cosh/p/2 . 

As all elements of modular group matrices are integers the trace is also an 
integer 

|TrM| = n . 

In Sect. ^ it was demonstrated that the mean multiplicity of periodic orbit 
length for the modular group is 

J/2 

(5(0) = 2— . 

Denote by n the trace of a given conjugacy class and by g(n) the number of 
distinct conjugacy classes corresponding to trace n. As n goes as e^^^ when 
n ^ oo one concludes that 

(gin)) ^ . (80) 

inn 

According to the Selberg trace formula the density of eigenvalues for the 
modular surface d{E) = d{E) + d^"^^' [E) where the oscillating part of the 
density is represented by the following formal sum 

S°^''\e) = Ay 3(")— cos(2A:lnn) . 

n 

From (|80|) it follows that mean value of g{n) Inn/n is one. Therefore we define 

Inn 

a[n) = g{n) , 

n 

so 

d^°"'\E) = ^ Va(n)cos(2fclnn) 
Trfc ^-^ 

n 

and (a(n)) = 1. 

As it was done in Sect. ^ one gets 

i?2(ei,e2) =^P + i?2(ei,e2) 



where 
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- \ n( 71 

(27rfc)2 



ni,n2 
2i(fci Inni — fe2 lnn2) 



and 



As was discussed in Sect. ^ due to the energy average the first term wih be 
washed out and the second one gives contributions only when 

ri2 = ni + r with r <^ni . 

Finally R^i^i, £2) — R2{f) where e = ei — £2 and 

— 1 \ ^ \ ^ / N / \ f e^-f^''" ■ lii"- 

^2 (e) = 2^ 2^ a(n)a(n + r) exp [ -2i— + le 



Let assume that the following mean value exists 

N 



(r) — lim — a(n)a(n + r) 



7, . , 

The dominant contribution to the two-point correlation function corresponds 
to replace the product OL(n)a{n + r) by its mean value 7(r) 



s=w-ji./;*.f:.He-»'-'"»p(i.!^) 

"''^ r— — 00 ^ ^ 



(81) 



where we have used a continuum approximation for n starting formally from 
a certain fixed rip ^ 1 , since only large values of n make a significant contri- 
bution. 

Define a (real) function /(x) as follows 

00 

!{x) = 7Me-^" . (82) 

r— — 00 

This function has the meaning of the Fourier transform of the two-point cor- 
relation function for multiplicities of the modular group. 

After the changing variable n e^^ in H81|l one gets that the two-point 
correlation function for the modular group is expressed through f{x) as follows 

— 1 f°° 

R2{e)^^rr e'="/(2fce-'^'") coseu dw (83) 
and the two-point correlation form factor is 
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K{t) = / i?2(e)e-*de = e'^' f{2ke-''') . (84) 

Therefore all non-trivial information is contained in functions 7(r) or /(x). 

The simplest diagonal approximation is to assume that the a{n) are essen- 
tially uncorrelated, that is, 7(r) is zero for r ^ 0. This gives for f{x) a constant 
value which leads to an exponential growth of K{t) as in (|79|l . But from gen- 
eral considerations K(t) obtained from a discrete spectrum has to saturate to 
a constant value for i ^ oo, consequently, the diagonal approximation cannot 
be correct for large t. 



4.2 Two-Point Correlation Fiinction of Multiplicities 

The purpose of this Section is to calculate the two-point correlation function 
of modular group multiplicities, j{r), whose Fourier harmonics according to 
H84(l determines the two-point correlation form factor. 

The calculation will be done by a generalization of the Hardy-Littlewood 
method for prime pairs discussed in Sect. 12. l1 As for primes one has to perform 
the three following steps. 

The first step 

Define the mean value of a(n) when n runs over integers of the form mq + r 
for fixed q and r < g in the following way 

1 

a{q;r) — lim — a(rnq + r) . 

m=0 



Since {a{n)) — 1 



9-1 

"^a(g;r) = q 



r=0 



Let Mq be the set of 2 x 2 matrices with entries being integers modulo q and 
having determinant equals one modulo q. These matrices form a group under 
multiplication modulo q which is sometimes called the modulary group. 

Define also Mg^r to be the set of elements of Alq with trace equal to r 
modulo q. One can prove that 



a{q;r) 



q\Mq.r\ 



where |M| is the number of elements of a set M. 

The intuitive meaning of this result is the following: g{n) is the number 
of conjugacy classes of modular matrices of trace n. To each modular matrix, 
one can associate an element of Mq in a unique way simply by taking the 
entries of the matrix modulo q. If n is equal to r modulo q, then all these 
matrices will belong to Mq^r- If we therefore assume that the matrices of the 
modular group cover the set Mg in some sense uniformly, the result follows. 
More careful treatment has been performed in j28| . 
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Example 

Let us consider q = 2. Integers modulo 2 are and 1. The group M2 consists 
of the following matrices 

10\ /10\ /01\ /01\ /II 

OlJ '\OlJ J '\10 J J '\10 

The dimension of the group M2, i.e. the total number of matrices, IM2I = 
6. Besides these matrices there are four matrices with zero trace (mod 2), 
i.e. |M2,o| = 4, and two matrices with trace equals 1 (mod 2), |M2,i| = 2. 
Therefore 

2-44 2-22 

The second step 

Define as in the Hardy-Littlewood method the following function 



n=0 



Since {a{n)) = 1 the convergence radius of this series is equal to one. 
The importance of this function follows from the integral 



J^(e-") = e™ / (e-"+'^) # (e"""'^) e"^''^ = V a{n)a{n + r)e-2"" 

Jo 27r 



whose right-hand side by a Tauberian theorem is connected with the two-point 
correlation function of multiplicities, j{r). 

The essence of the Hardy Littlewood approach is the investigation of the 
function ^{z) when z = exp(— m -|- ie -|- 2'K\p/q) with u ^ and e ^ 0, where 
p and q are co-prime integers. The main step is then to write n in the form 
mq + r with r lying between and q — I and prove that in the expression 
for ^{z) the dominant contribution as u and e go to zero will be given by the 
mean value of a(mg -|- r), that is, one may substitute it by a{q; r). 

Accepting this, one has that as u — > and e — > 



q—l 00 



'P (exp {-u + 2TTip/q + ie)) = XI "(^9 + r)e-("-'<^)('"«+'')e2'''''P/« 

r=0 m=0 

= Va(g;r)e2'^'3"'/«- / dne-^''-^^)" = 



r=0 

u — ie 



where 
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9-1 



■ P 

' r 



f3{p,q)=q ^^a((7;r)exp f 27ri 

Hence ^(z) has a pole singularity at all rational points on the unit circle. 
The third step 

Divide the unit circle in intervals Ip,q centered around exp(27rip/(7), where 
p and q are co-prime integers with p < g. If one neglects all terms in each 
interval except the pole term and extends the integration over e to the whole 
hne, one gets 



J.(e-")=e™ I 



(P,9) = l' 



de\f3{p,q)\^ 

27r u2 + ~ 



^ Y mp,q)fexp(2m^r 

2u ^-^ V q 



Finally one obtains that 

= Y^ l/3(p,9)l^exp f 27ri-r 
(p,«)=i ^ ^ 

The sum is performed over all and p co-prime to q with Q < p < q. 

This is the two-point correlation function of multiplicities of the periodic 
orbits for the modular group . All other quantities of interest can be obtained 
from it. In particular, the function ](x) l|82|) is given by 

/(x) = 27r Y l/3b,'z)P<5(x-27r^ 

(P,9) = l ^ ^ 

where the summation is done over all p and q co-prime, without the restriction 
p<q. 

According to (|83|l and (|84|l the knowledge of ](x) determines immedi- 
ately the two-point correction function the form factor of modular domain 
eigenvalues. 

4.3 Explicit Formulas 

Let us define the so-called Kloosterman sums 

5(n, m;c)= "Y^ 6xp | — -(nd + mc?^^) 



(d,c) = l 



where the summation is taken over all d < c co-prime with c and d -'^ is an 
integer modulo c which obeys d~^d~ 1 (mod c). 
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One can show (see that q) can be expressed through these sums 
in the following way 

where lo are the prime divisors of q. 
The function 7(r) can be written as 

oo 

where Ar{q) is given by 



One can prove that Ar{q) is multiplicative function of q, i.e. Ar{nin2) = 
Ar{ni)Ar{n2) provided {nin2) = 1, therefore one needs to know only its 
values on powers of primes and j{r) can be rewritten as the infinite product 
over all prime numbers 



k=l 



To present a closed expression for Ar{q) let us introduce the standard defini- 
tion of the Legendre symbol 



1 , if a = a;^ (mod q) has a solutiona ^ (mod q) 
, if a = (mod q) 
'^^ 1^1' otherwise 

The meaning of this symbol is perhaps best understood by saying that the 
number of distinct solutions of the equation x'^ = a (mod is 1 + (a/q). 

A fairly tedious evaluation of Ar{q) (see for details) gives the following 
formulas. 

Let q = p^ where p is an odd prime. Then for n = 1 one has 



Arip) = 



(p2 „ 1)2 



•^g/ (x^-4)((x + r)^-. | ^ 



x=0 



For n > 2 we have, letting t be an arbitrary non-zero number modulo p 

Aripn 



2(1 - 1/p) , r = (mod p") 

1 I -2/p , r EE (mod p") 

p2„(i _p-2) ^ e(n,p)(l - 1/p) , r EE ±4 (mod p") 

~e{n,p)/p , r = ±A + tp^^^^ (mod p") 
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where e{n,p) takes the value —1 if n is odd and p is of the form 4fc + 3 and is 
equal to 1 in all other cases. For p = 2, we list down individual cases for low 
powers and eventually state a general rule 



Ar{2) = 

A. (4) = 
A. (8) ^ 

Ar{lQ) = 





1 , 


r = 




-1 , 


r = 1 




' 1 


, r = 


18 ^ 




, r = 2 


, 







(mod 2) 
(mod 2) 

(mod 4) 
(mod 4) 



1 



9-16 

A^(32) = , 

and finally, for the general case rt > 6 







(mod 16) 
(mod 16) 



1 



9-22" 



r EE 

— nn— 1 



2^' 

±(4 + 2"-2) 

±(4 + 2"-2 + 2"-i) 



(mod 2") 
(mod 2") 
(mod 2") 
(mod 2") 



All terms not explicitly shown equal zero. In |50| these formulas were proved 
by a different method. 



4.4 Two-Point Form Factor 

These formulas give the explicit expression for the two-point correlation form 
factor 

1 



2n^k 



E 

(P,9) = l 



P 



Sit - tp^q) 



where 



1, kq 

fc np 



In the limit fc — > cx) and t fixed, the dominant contribution comes from terms 
with p/q <^ 1. Smoothing over such values one can show (see that in this 
limit K{t) tends to the constant Poisson value 

A 



Kit) = 



47r 



where A = 7r/3 is the area of the fundamental region of the modular group. 
For small t (of the order of the Ehrenfest time \nk/k) Kit) has number- 
theoretical oscillations due to cumulative contributions of degenerate periodic 
orbits. For very small values of t (of the order of 1/fc) the two-point form 
factor has S peaks connected with short periodic orbits. 

Though the modular group is by no means a generic system, it is the first 
ergodic dynamical system for which it was possible to compute explicitly the 
distribution of the energy levels. 
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5 Hecke operators 

Arithmetic groups have many interesting properties. In particular, for aU 
arithmetic groups it is possible to construct an infinite number of mutually 
commuting operators which commute also with the Laplace-Beltrami oper- 
ator. These operators are of pure arithmetic origin and are called the Hecke 
operators |37]) Ell- 
in a certain sense these operators permit to 'understand' why arithmetic 
systems have the Poisson statistics typical only for integrable systems. The 
point is that integrable systems are systems with sufficiently large number 
of independent commuting operators and Hecke operators may be viewed as 
a manifestation of a kind of arithmetic integrability of arithmetic systems 
which does the Poisson statistics for these models natural 123 . Unfortunately, 
precise relations along this line seem to be impossible. 

Let us consider informally the construction of Hecke operators for the 
modular group. Choose two matrices A and B from the modular group with 
the same trace. As they have the same trace and determinant, they have the 
same eigenvalues and there exists a matrix 7 such that 

7^7-1 ^ B or -/A^ B-f and det(7) ^ . (85) 

If A and B are not congugated in the modular group, 7 ^ PSL(2, 7Z). But 
matrix 7 can be chosen as a matrix with integer elements but with the deter- 
minant 1. 

Example. 

Consider the following simple matrices 

--{ID- 

General form of matrices 7 which obey (|85|l is 

_ /2a + 2/3 a \ 

with arbitrary a and (3. 

It is clear that there exists no 7 g PSL(2, TZi) but choosing different integer 
values of a and /3 one can construct an infinite number of integer matrices 
with determinant 7^ 1 which obeys (|85|l . For example, 

7- (det = -1), 7= (I J) (det = 3)... . 

These considerations demonstrate that when dealing with the modular group 
it is quite natural to consider matrices with integer elements but with the 
determinant different from one 
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Mp ~ "I ^ ^ a,b,c,d integers, ad—bc~p 

Matrices Mp with p ^ I do not form a group because their product has not 
the same form. 

A matrix mp € Mp can uniquely be represented in the form 

TTlp = /ittp (86) 

where fi G PSL(2, TZi) and ap is one of matrices from the foUowing finite set 
ap = I a,b,d integers , ad = p , d> , < b < d — 1^ . (87) 

Instead of proving this fact let us transform a simple matrix 



TO3 



4 1 
1 1 



to the form (|86|l . General proof (see e.g. ^54J follows the same steps. 
First, it is necessary to find a matrix 



such that (i) det n' = 1 and (ii) 



^i/X2\/4l\/a5 
^i3^^4J\llJ \0d 

The condition of the zero of low-left element gives the equation 4/i3 + /i4 = 
and because /i4 and /is are coprime they can be chosen as follows: (14—4, and 
/i3 = —1. Unit determinant condition gives fix — k and /i2 = 1 — 4fc with an 
arbitrary integer k. Finally, b — 1 — and the smallest positive b modulo 3 
corresponds to = 0. Hence, our matrix ma has the following representation 

An important property of the set Mp is that when one multiplies a matrix 
from this set by a matrix from the modular group the resulting matrix also 
belongs to Mp 

Mpg ^ Mp for all g £ PSL(2, ffi) . 
Let ^{z) be an automorphic function of the modular group, i.e. 

^{gz) = ^{z) , for all g e PSL(2, TL) . 

Then it is easy to see that the function 
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where the summation is performed over all ad — p, d > 0^0 < b < d ~ 1 will 
also be an automorphic function for the modular group. This is a consequence 
of the fact that in the right-hand side of this expression there is effectively 
the summation over all matrices from Mp. As Mp goes not change after mul- 
tiplication by a modular group matrix ^'{z) is an automorphic function for 
the modular group. {TpW)(z) is a kind of symmetrization of over images 
of z by all elements of Mp and the operators Tp are called Hecke operators. 

These operators form a commutative algebra with the following product 
(see e.g. jHH) 

TnTm — E Tnm/d^ (88) 
d\ ln,m,) 

where the summation is done over all divisors of the greatest common divisor 
of TO and n. The most important case corresponds to Hecke operators with 
prime indices because all the others can be obtained from H88|l. 
When p is a prime number 



{Tp^)[z)^^ 
VP 



W{pz) + E If' ^ ^ ^ 

0<j<p 



Since Hecke operators involve only fractional transformations all them com- 
mute with the Laplace-Beltrami operator. Consequently, if !f'(x, y) is an eigen- 
function of the Laplace-Beltrami operator, then {Tp^){x^y) will also be an 
eigenfunction with the same eigenvalue. If there is no spectral degeneracy 
(which strongly suggested by numerics) every eigenfunction of the Laplace- 
Beltrami operator is in the same time an eigenfunction of all Hecke operators 

{Tp^^){x,y) ^ \p{n)^n{x,y) . (89) 

It is known (see e.g. [SJ) that eigenfunctions of the Laplace-Beltrami operator 
for the modular group have the following Fourier expansion 

oo 

W^{x,y)^y^'^ E Cp(n)if,„_i/2(2^ra)e2"f- 

p— — oo 

where the eigenvalue of the Laplace-Beltrami operator £"„ = s„(s„ — 1) and 
K^{x) is the Hankel function. 

One has z = x + iy and {az + b)/d ~ {ax + b)/d + iay/d, therefore 



^ a.h.d p 



where the first summation is performed over all a, 6, d as in (|87|1 . 
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The summation over b gives zero if d is not divide p. Otherwise 

d\p,d\7n 

Let k — m/d and u = pm /(P . Then p — mu / and 

M fc|(m,u) 

If T,n^n — ^m{n)^n then by comparing the first Fourier coefficient one gets 

Cm{n) = A™(n)ci . 

Assuming ci ^ Q and using a convenient normalization ci — \ one concludes 
that eigenvalues of the Hecke operators coincide with the Fourier coefficients. 

We note also that similarly to the construction of the Selberg trace for- 
mula one can build the trace formulas for Hecke operators (see e.g. [23 and 
references therein). Such trace formula schematically has the form (cf. H29|l ) 

+ smooth, parabolic, and elliptic terms . 

Here h{k) is a test function like in Sect. 12.81 and g{l) is its Fourier transform. 
In the left-hand side the summation is performed over all eigenvalues i?„ — 
+ oi the Laplace-Beltrami operator and Ap(n) is the eigenvalue of the 
Hecke operator Tp (|89|l applied to the eigenfunction of the Laplace-Beltrami 
operator with eigenvalue En- In the right-hand side the summation is done 
over all 'hyperbolic' matrices from Mp with Tr rrip ^ p+1. Lp is the 'length' 
associated with matrix mp 

2cosh(ip/2) = I Tr mp|/VP 

and lp is the minimal length of modular group matrices commuting with mp. 

6 The Jacquet— Langlands Correspondence 

Another curious fact about arithmetic groups is the Jacquet-Langlands cor- 
respondence (see 40 ) which claims that for a arithmetic group derived from 
a quaternion group overQ (with a finite fundamental domain) one can find a 
subgroup of the modular group (with infinite fundamental domain) in such a 
way that amongst all automorphic eigenvalues of the Laplace-Beltrami opera- 
tor for this modular group subgroup one can find all eigenvalues of a compact 
arithmetic group. 
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The simplest arithmetic group F derived from quaternion algebra over Q 
with division is (see Sect. I2.2|l 



r 



f fci + k2\/a + ki^/a 
\h{k^ - k^y/a) ki - k2\fa 



where 6 is a prime number, a is an integer such that the equation = 
a ( mod 6) has no integer solution (e.g. a = 3, 6 = 5), and integers ki are such 
that 

det(7) = k\- ak\ - hk\ + ahk\ = 1 . 
Denote z — x -V iy^ t ~ u ^ iv (y,v > 0) and define for all rij 

a = ni + n2\fa , j3 = + rii^/a , 

7 = 6(^3 — ni^/a) , S ~ ni — n^yfa . 
Fix an arbitrary zo and compute the following kernel 

Ht,z)= expK{T,z) 

Uj — — OO 

where 

K{r, z) = -.Imr '""° + f-f^'° + ^^'' + 2.rf {aS - Pj) . 

Im zlmzQ 

Here z is the complex conjugate of z. 

Let ijjniz) be an eigenfunction of the Laplace-Beltrami operator automor- 
phic with respect to the quaternion group F. It means that 

• il)n{Mz) = -tpn^z) for ah M e F. 
Then the function 

dxdy 



^{t) = IfjiT / <Z'(r,z)V'„(2;) 
Jv 



where the integral is taken over the fundamental domain of the group F is an 
eigenfunction the Laplace-Beltrami operator with the same eigenvalue £"„ but 
automorphic with respect to the congruence subgroup of the modular group 
/o(4afe) where 

n(A^)= (X") eSL(2,K) 
with an additional condition that 

fc = ( mod N) . 
Direct (but tedious) proof of this statement can be found in j 4()| . 
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7 Non-arithmetic Triangles 

In the precedent Section we have seen that arithmetic systems have the Pois- 
son spectral statistics. But what is about non-arithmetic models? 

Let us consider, as example, the so-called Hecke triangles which are hy- 
perbolic triangles with angles (0, 7r/2, 7r/n). All of them tessellate the upper 
half-plane and are fundamental domains of the discrete groups generated by 
reflections across their sides. The modular billiard is a part of them corre- 
sponding to n = 3. Similar to it they all have an infinite cusp. 

According to Table the Hecke triangles are arithmetic only for n — 
3,4, 6, cx). All these arithmetic triangles have an exponential degeneracies of 
periodic orbit lengths which leads to the Poisson-like statistics of energy levels. 

The simplest non-arithmetic Hecke triangle is the one with n — 5. At 
Fig-Elwe present the results of numerical calculations of the nearest- neighbor 
distribution for 6000 first energy levels for this triangle with the Dirichlet 
boundary conditions. For others Hecke triangles one gets similar pictures. 
It is clearly seen that numerics agrees very well with the predictions of the 
Gaussian Orthogonal ensembles of random matrices as it should be for generic 
chaotic models. 




12 3 4 



Fig. 19. The nearest-neighbor distribution of 6000 energy levels for the non- 
arithmetic Hecke triangular billiard with n = 5. The solid line - the GOE prediction. 
Dotted line - the Poisson result. 

But what are the multiplicities of periodic orbit lengths for non-arithmetic 
Hecke triangles? As these model are not-arithmetic, one would expect that 
their length multiplicities should be equal to two as for generic time-reversal 
invariant systems. Nevertheless numerical calculations (see for details) 
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demonstrated that this is not always the case. At Fig. 1201 we present the 




5 10 15 20 

/ 



Fig. 20. Mean length multiplicities of periodic orbits for the Hecke triangles with 
(from top to bottom) n — 12, n = 5, n = 8, and n = 10. White lines are numerical 
fits TO . 

numerically computed mean length multiplicities for the Hecke triangles with 
n = 5, 8, 10, 12 for lengths / < 20. White lines indicate a two-parameter fit 
to these data in the form g{l) = a„e''"' 

n = 5 : g{l) w 1.235e "'*' , 

n = 8 : g{l) w LOgSe "-*' , (90) 
n = 10 : g{l) w 1.143e °^^' , 
n = 12 : g{l) « .986e-^5"' . 

These expressions fit numerical data in the given interval of lengths quite well 
and indicate that for, at least, certain Hecke triangles mean length multiplic- 
ity increases exponentially. We stress that (|90(l are only the best least-square 
numerical fits and no attempts were made to determine the accuracy of coef- 
ficients. 

The discussion of the origin of such unexpected multiplicities for non- 
arithmetic triangles is beyond the scope of these lectures (on this subject 
see However it is of interest to understand why exponentially large 

multiplicities of periodic orbit lengths do not contradict the observed GOE 
behaviour of spectral statistics (cf . Fig. I19|) . 

Assume that a system has an exponentially large number of periodic orbits 
with the same length I increasing as 
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with A < 1/2. 

Let us repeat the arguments of Sect. 11.21 for this case with exact degen- 
eracicies. In Sect. II .21 it was demonstrated that periodic orbits with different 
lengths can be treated in the diagonal approximation if 

k 

Ipi ~ ^ (91) 

where k is the momentum and AE is the width of the energy average inherent 
in the definition of correlation functions of dynamical systems. 
As the density of orbits with different lengths is 

e' 

Pdiff. lengths ~ ~(T\] ^ C^"^""^^' 

it follows that the inequality H91|) is valid till maximal length 

/„.^^ln^^^ln.. (92) 

Notice that due to assumed large multiplicity Im is different from (|^ . 

From H78|) it follows that the two-point correlation form factor in the di- 
agonal approximation up to numerical factor is 

i^(t)^y|A(0|W 

where t = l/2k and A{1) ~ y/^/k. Comb ining all terms together one obtains 
that during the maximal time of applicability of the diagonal approximation 
tm = Im/'^k with l„i from (|92ll the form factor increases till 

k 

Hence, if A = 1/2 as for arithmetic systems the two-point correlation form 
factor during the time of validity of the diagonal approximation increases 
till a constant value of the order of 1. But if A < 1/2 the form factor for 
the time of validity can reach only a value of the order of k^'^ with v — 
(1 — 2A)/(1 — A) > 0. As fc ^ oo this value tends to zero and no apparent 
contradiction with standard random matrix ensembles can be derived within 
the diagonal approximation. 



8 Summary 



Arithmetic groups are a special sub-class of discrete groups characterized by 
the existence of a representation by matrices with integer elements. A readable 
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mathematical review of such groups is given in ^42. . There are two types of 
arithmetic groups. The first includes groups commensurable with the modular 
group and having non-compact fundamental domains with infinite cusps. The 
second type of compact arithmetic groups combines groups commensurable 
with groups derived from quaternion algebras with division. These groups 
have finite fundamental domains. 

From classical viewpoint the free motion on surfaces generated by arith- 
metic groupes is as chaotic as for any hyperbolic surfaces. But quantum me- 
chanics on these arithmetic surfaces is very special. In particular, spectral 
statistics of the Laplace-Beltrami operator automorphic with respect to arith- 
metic group is described by the Poisson statistics typical for integrable systems 
and not by the random matrix statistics typical for chaotic models. 

The origin of this peculiarity can be traced to the existence in arithmetic 
systems of a very large number of periodic orbits with exactly the same length. 
For all arithmetic groups the mean multiplicity of periodic orbits with length 
/ behaves like e'^^/L This has to be compared with the total density of peri- 
odic orbits which for all discrete groups is /I. It is the cumulative effect of 
the interference of many periodic orbits with the same length which changes 
drastically the spectral statistics. 

In the diagonal approximation the two-point correlation form factor K{t) 
for arithmetic systems at small t increases exponentially like e'^'/A: and dur- 
ing the Ehrenfest time (which is the limit of applicability of the diagonal 
approximation) reaches a constant value. 

More detailed information can be obtained for the modular group where it 
is possible to compute the two-point correlation form factor analytically. The 
final answer is 

2 

where 

1, kq 

'tp,q = 7 111 — 
K irp 

and (3(p, q) is a number-theoretical function given in Sect. 14.31 

This formula means that the two-point form factor for the modular group 
is a sum over (5-functions at special points tp^q situated in a vicinity of the 
Ehrenfest time. The set of ^-functions is dense but the largest peaks corre- 
spond to the smallest ratios p/q. Nevertheless, small peaks with p/q <C 1 are 
much more numerous and integrally they dominate. In the limit t fixed and 
k —> oo K{t) —f d thus confirming the Poisson nature of the spectral statistics 
of the modular group. 

Arithmetic groups have many interesting properties. Hecke operators and 
the Jacquet-Langlands correspondence are the most remarkable. 



K{t) = 



27r2fc 



E 
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